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of  the  Requirements  for  the  Degree  of  Doctor  of  Philosophy 

A STUDY  OF  THE  GLOBAL  EVOLUTION 
OF  MICROSTRUCTURES 

By 

Arun  Mahadeo  Gokhale 
March,  1977 

Chairman:  R.  T.  DeHoff 

Major  Department:  Materials  Science  and  Engineering 

Mechanical  and  physical  properties  of  crystalline 
materials  depend  on  the  geometrical  properties  of  their 
microstructures.  Thus,  an  understanding  of  the  processes 
by  which  microstructures  evolve  is  of  great  significance. 

In  many  cases,  theoretical  and  experimental  studies  on  the 
physics  of  microstructural  evolution  are  done  on  a "local” 
basis.  On  the  other  hand,  the  experimental  techniques  of 
Stereology  usually  yield  the  "global"  or  average  geometri- 
cal properties  of  microstructures.  Thus,  there  is  a need 
to  bridge  the  gap  between  the  local  and  the  global  de- 
scription of  the  evolution  of  microstructures.  This  is 
the  major  objective  of  this  dissertation. 

An  approach  has  been  developed  to  evaluate  the  global 
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geometrical  properties  of  an  evolving  microstructure  from 
a given  description  of  the  local  variables.  Effects  of 
various  growth  geometries  on  the  kinematics  of  phase 
transformations  have  been  explored  in  detail.  A pheno- 
menological approach  is  developed  to  evaluate  the  impinge- 
ment equations  for  two  and  one  dimensional  growth  geo- 
metries. This  phenomenological  approach  is  valid  for 
both  random  and  non-random  impingements.  A fairly 
general  technique  is  developed  to  deduce  the  nucleation 
rate  and  growth  rate  from  the  time  variation  of  the 
global  properties. 

The  practical  feasibility  of  the  approach  is 
demonstrated  by  applying  it  to  the  experimental  data  on 
recrystallization  and  austenitization  processes.  The 
analysis  is  used  to  deduce  the  nucleation  rate  and  growth 
rate  of  austenite,  during  the  austenitization  of  a 
pearlitic  structure.  It  is  concluded  that  the  nucleation 
rate  and  growth  rate  of  austenite  remain  constant  at  a 
given  temperature  during  austenitization  of  pearlitic 
structures . 

The  kinetics  of  austenitization  of  some  "spherodized" 
cementite  structures  is  studied  in  detail.  It  is  con- 
cluded that  the  "spherodized"  cementite  is  present  in  the 
form  of  an  incomplete  grain  edge  network,  which  dissolves 
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in  a one  dimensional  manner  during  the  austenitization 
process.  The  austenite  grains  grow  along  the  ferrite 
grain  boundaries  in  a two  dimensional  manner.  The  present 
approach  is  successfully  applied  to  explain  the  experi- 
mentally observed  behavior  of  nine  independent  stereolo- 
gical  properties. 


X 


CHAPTER  I 

PROBLEM  FORMULATION 


The  properties  of  materials  are  to  a large  extent 
determined  by  their  ''structure.”  In  crystalline  mate- 
rials , atoms  or  molecules  of  the  constituents  are  arranged 
in  a periodic  array  in  space.  This  gives  rise  to  crystal 
structure.  The  same  material  in  two  crystalline  forms  may 
have  very  different  properties.  On  a coarser  level,  the 
structure  is  classified  in  terms  of  heterogeneity  of  its 
individual  crystals  or  grains.  Heterogeneity  arises  be- 
cause more  than  one  crystal  may  nucleate  and  grow  when  a 
material  crystallizes  or  transforms  from  one  crystalline 
form  to  another.  This  gives  rise  to  microstructure. 
Mechanical  and  physical  properties  of  materials  depend  on 
the  geometrical  properties  of  their  microstructures.  For 
example,  yield  stress  varies  linearly  with  the  square  root 
of  the  mean  grain  intercept  [1].  It  may  be  said  that  con- 
trol of  microstructure  would  lead  to  the  control  or  opti- 
mization of  some  property  or  a set  of  properties.  In 
order  to  control  microstructures,  it  is  necessary  to  con- 
trol the  processes  by  which  microstructural  changes  occur. 
Thus  the  understanding  of  the  processes  by  which  micro- 
structures evolve  is  of  great  significance. 
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In  many  cases,  theoretical  and  experimental  studies 
on  the  physics  of  microstructural  evolution  are  done  on  a 
"local"  basis.  For  example,  diffusion  controlled  growth 
of  a precipitate  is  described  in  terms  of  the  local  inter- 
face velocity  or  local  growth  rate.  These  theoretical 
models  describe  how  a microstructure  evolves  locally  [2- 
11].  On  the  other  hand,  quantitative  structure-property 
correlations  are  usually  obtained  through  the  "global"  or 
average  properties  of  microstructures.  Furthermore,  the 
experimental  techniques  of  quantitative  microscopy  usually 
yield  only  the  global  or  average  geometrical  properties  of 
microstructure.  Thus,  there  is  a need  to  bridge  the  gap 
between  the  local  and  the  global  description  of  the  evolu- 
tion of  microstructures.  This  is  the  major  objective  of 
this  thesis. 

The  "local"  variables  and  the  kinetic  variables  which 
affect  the  "global"  geometrical  properties  of  a micro- 
structure, during  and  at  the  end  of  a phase  transformation 
or  a redistribution  process,  are  as  follows: 

(1)  Nucleation  rate 

(2)  Growth  or  shrinkage  rate 

(3)  Particle  shape 

(4)  Growth  or  shrinkage  geometry 

In  addition  to  these  local  and  kinetic  variables , the 
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global  microstructural  evolution  is  determined  by  the  fol- 
lowing variables: 

(5)  Starting  microstructure 

(6)  Distribution  of  nuclei  in  space 

The  answers  to  the  following  questions  provide  the  solu- 
tion to  the  problem. 

(A)  If  a quantitative  description  of  these  six  variables 
is  given,  how  are  the  geometrical  global  properties  of  a 
microstructure  during  its  evolution  evaluated? 

(B)  Given  some  global  geometrical  properties  of  micro- 
structure at  different  times  during  a phase  transformation 
or  a redistribution  process,  how  to  evaluate  the  local 
variables  ? 

An  attempt  will  be  made  to  answer  these  two  questions 
in  the  third  chapter  of  this  thesis.  However,  some  phys- 
ically acceptable  simplifying  assumptions  will  be  made  to 
solve  the  problem.  This,  then,  leads  to  the  following 
question : 

(C)  Is  such  an  approach  practically  feasible? 

In  order  to  demonstrate  the  feasibility  of  the 
approach,  it  will  be  applied  to  the  experimental  data  on 
the  recrystallization  and  austenitization  processes.  The 
fourth  chapter  deals  with  the  interpretation  of  the  ex- 
perimental data  on  the  recrystallization  in  some  Fe-Si 
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alloys.  The  fifth  chapter  deals  with  the  experimental 
measurements  of  the  global  properties  of  austenite  during 
the  austenitization  of  pearlite  and  the  evaluation  of  the 
local  variables  from  the  time  variation  of  the  global 
properties  of  austenite.  The  sixth  chapter  is  concerned 
with  a detailed  quantitative  microscopic  investigation  of 
the  austenitization  of  spherodized  cementite  structures. 
The  present  approach  will  be  successfully  applied  to 
interpret  the  experimental  data. 


CHAPTER  II 

STEREOLOGICAL  METHODS 


2 . 1 Introduction : 

In  virtually  every  branch  of  science,  there  exists  an 
interest  in  the  quantitative  characterization  of  struc- 
tures that  appear  in  nature.  The  method  of  observation 
of  such  structures  always  involves  some  biased  sampling 
technique.  For  example,  observation  of  a microstructure 
with  an  optical  microscope  presents  a two  dimensional 
section  through  the  three  dimensional  structure  contained 
in  the  sample.  Measurements  made  upon  such  a two  dimen- 
sional section  are  related  to  the  geometry  of  the  three 
dimensional  structure  which  they  represent.  The  process 
of  quantitative  evaluation  of  some  properties  of  the  three 
dimensional  microstructure  is  called  "Quantitative  micro- 
scopy" or  "Stereology . " Apart  from  Metallurgy,  quanti- 
tative microscopic  techniques  are  used  in  a variety  of 
fields  such  as  Geology  [12,13],  Biological  Sciences  [14, 
15],  Ceramics  [16],  Concrete  technology  [17],  Wood  tech- 
nology  [18],  etc.  What  follows  is  primarily  a catalog  of 
the  basic  relationships  of  quantitative  microstructural 
analysis . 
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2.2  Measurement  of  Global  Properties : 

The  global  properties  of  a microstructure  are 
classified  into  two  categories,  namely,  metric  and  topolo- 
gical properties.  The  metric  properties  (for  example, 
volume,  surface  area,  etc.)  can  be  evaluated  from  the 
measurements  made  on  a two  dimensional  section,  whereas 
the  topological  properties  can  not  be  evaluated  from  such 
measurements  without  making  any  simplifying  assumptions. 

A serial  sectioning  of  the  structure  is  necessary  for 
evaluating  the  topological  properties  of  an  opaque  body. 
The  metric  properties  describe  the  numerical  extent  of  a 
particular  feature  in  the  structure.  Each  three,  two  and 
one  dimensional  feature  in  the  structure  has  associated 
with  it,  respectively,  volume,  area  and  length,  which  can 
be  measured  unambiguously  from  the  plane  of  polish.  The 
procedure  for  such  measurements  will  now  be  described 
briefly. 

2.2.1  Measurement  of  Volume  Fraction 

The  volume  fraction  of  a particular  phase  in  the 
structure  can  be  measured  by  using  areal  analysis  [19-24], 
linear  analysis  [20-22,24-26]  or  point  counting  [20,27]. 

The  basic  principles  involved  in  these  three  methods  are 
as  follows. 


7 


Areal  analysis:  The  expected  value  of  the  area 

fraction  of  a particular  phase  measured  on  the  plane  of 
polish  is  identically  equal  to  its  volume  fraction  in  the 
three  dimensional  structure. 

Lineal  analysis:  Suppose  a number  of  test  lines  are 

superimposed  on  the  plane  of  polish;  then  the  expected 
value  of  the  fraction  of  length  of  the  test  lines  con- 
tained in  a particular  phase  is  identically  equal  to  its 
volume  fraction. 

Point  counting : If  a point  is  placed  at  random  on  a 

random  plane  of  polish,  then  the  probability  of  the  point 
falling  within  a particular  phase  is  equal  to  the  volume 
fraction  of  that  phase. 

A detailed  analysis  of  these  three  techniques  is 
given  by  DeHoff  and  Rhines  [20]  and  Underwood  [26]  and 
will  not  be  repeated  here. 

2.2.2  Measurement  of  Surface  Area  per  Unit  Volume 

A statistically  exact  relationship  for  the  evaluation 
of  interfacial  area  per  unit  volume,  Sy,  was  first  de- 
veloped by  Saltykov  [22]  in  1945  and  later  by  Smith  and 
Guttman  [28].  Alternative  derivations  are  given  by 
Underwood  [26]  and  DeHoff  and  Rhines  [20].  Suppose  a 
number  of  test  lines  are  superimposed  on  the  plane  of 
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polish.  Let  'L'  be  the  total  length  of  these  test  lines 
and  'N'  be  the  number  of  intersections  of  these  test  lines 
with  a particular  type  of  interface  (e  g.  grain  boundaries, 
P^^ticle-matr ix  interface,  etc.),  then  the  surface  area 
per  unit  volume,  S^,  of  that  interface  is  given  by 

Sy  = 2*Expected  value{Nj^}  (2.1) 

where  = N/L 

If  directional  anisotropy  exists  in  the  structure,  then 
different  orientations  of  the  test  lines  will  give  dif- 
ferent values  of  Sy.  This  concept  was  used  by 
Hilliard  [29,30]  in  his  treatment  of  oriented  structures. 


2.2.3  Total  Curvature 

Consider  a small  surface  element,  dS,  of  an  arbi- 
trarily curved  surface.  The  local  mean  curvature,  H,  of 
this  surface  element  is  given  by  [31] 

H = ^[K^  + K2]  (2.2) 


where  and  K2  are  the  principal  normal  curvatures . 

The  total  curvature.  My,  of  a given  type  of  interface 
per  unit  volume  of  the  structure  is  given  by 


M.  = //  H dS 


(2.3) 


where  the  integration  is  carried  out  over  all  the  surface 
elements,  of  the  interface  under  consideration,  in  a unit 
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volume  of  the  structure.  For  a system  of  discrete  convex 
particles,  the  total  curvature  is  proportional  to  the 
first  moment  of  the  size  distribution  function. 

DeHoff  [32]  has  shown  that  for  a tubule  network,  the  total 
curvature  is  proportional  to  the  length  of  the  network. 

The  total  curvature  of  surfaces  can  be  measured  quanti- 
tatively by  using  the  sampling  procedure  devised  by 
Rhines  [33]  called  the  area  tangent  count.  It  was  inde- 
pendently shown  by  DeHoff  [34]  and  Cahn  [35]  that  the 
total  curvature  of  an  interface  is  given  by 

= TT* Expected  value{(T.) 

= 2TT-Expected  } (2.4) 

^"^A^net  number  of  tangents  per  unit  area  formed 

by  a random  sweeping  test  line  with  the  traces  of  the 
surfaces  on  the  plane  of  polish.  Tangents  produced  with 
the  elements  of  traces  of  the  surface  that  are  convex  with 
respect  to  the  phase  of  interest  are  counted  as  positive 
and  those  produced  with  the  concave  elements  are  counted 
as  negative.  (^A^net  number  of  closed  area 

loops  per  unit  area  of  the  plane  of  polish.  Loops  which 
enclose  phases  other  than  the  phase  of  interest  are 
counted  as  negative. 
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2.2,4  Measurement  of  the  Geometrical  Properties  of 
Lineal  Features  ~ “ 

The  length  of  a lineal  feature  per  unit  volume  of  the 

structure,  Ly,  can  be  evaluated  by  using  the  following 

stereological  relationship  [22,28]: 


where  is  the  number  of  intersections  of  the  lineal 
feature  with  the  unit  area  of  plane  of  polish. 

DeHoff  [36]  has  shown  that  the  total  curvature  of  an 
edge  is  given  by  the  following  expression: 


The  integration  is  carried  out  over  the  length  of  the  edge 
in  the  unit  volume  of  the  structure  and  'x'  is  the  local 
dihedral  angle.  DeHoff  and  Gehl  [45]  have  shown  that  the 
average  dihedral  angle,  x>  is  given  by 


where  is  obtained  by  applying  the  area  tangent 

count  [33]  only  to  the  edges.  The  problem  of  the  measure- 
ment of  average  dihedral  angle  has  been  also  analyzed  by 
Harkar  and  Parker  [37] . 

2.3  Measurement  of  Particle  Size  Distribution: 


Ly  = 2 ‘Expected  value{P^} 


(2.5) 


f^V^edge  " ? / X 

■^TT 


V 


(2.6) 


(2.7) 


A particle  size  distribution  can  be  evaluated  from 
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the  measurements  made  on  the  plane  of  polish  if  some 
simplifying  assumptions  are  made  regarding  the  particle 
shape.  The  techniques  for  the  measurement  of  particle 
size  distribution  (henceforth  called  P.S.D.)  assume  a 
model  shape  for  all  the  particles,  and  require  a measure- 
ment of  the  distribution  of  some  geometrical  variable 
(e.g.  linear  intercepts,  section  sizes,  etc.)  on  the  plane 
of  polish.  This  distribution  function  is  then  converted 
into  the  three  dimensional  P.S.D.  by  using  the  principles 
of  geometrical  probabilities. 


2.3.1  Techniques  Based  on  the  Measurement  of  Section 
Size  Distribution 

Wicksell  [38],  Scheil  [39,40],  Schwartz  [41]  and 
Saltykov  [22]  have  developed  methods  for  evaluation  of 
P.S.D.  of  spherical  particles  from  the  distribution  of 
their  section  sizes.  Saltykov's  method  [22]  is  more  con- 
venient and  involves  less  computational  errors.  Excellent 
review  articles  on  the  subject  are  given  by  DeHoff  and 
Rhines  [20]  and  Underwood  [26]. 

All  these  methods  are  essentially  numerical  solutions 
to  the  following  integral  equation  [20] : 


“ n (R)dR 

f / T n — YJY 

^ r [R^-r^]^ 


(2.8) 


where  n. (r)  is  the  section  size  distribution  function  and 
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n^(R)  is  the 
tion  to  this 


three  dimensional  P.S.D.  An  analytical 
integral  equation  is  as  follows: 


n^(r)dr 


solu- 


(2.9) 


where 


N^^(R)  = / n^(R)dR  (2.9a) 

R 

The  direct  use  of  equation  (2.9)  requires  an  analytical 
form  for  the  section  size  distribution  and  furthermore  one 
has  to  evaluate  a derivative  if  the  three  dimensional 
P.S.D. , ny(R) , is  desired  instead  of  Ny^(R).  Due  to  these 
limitations,  the  numerical  solutions  [38-41,22]  to  equa- 
tion (2.8)  are  normally  used. 

In  the  general  case  of  ellipsoids  of  revolution,  with 
either  prolate  or  oblate  shape,  integral  equations  similar 
to  equation  (2.8)  can  be  set  up  [46].  DeHoff  [43]  has 
given  a solution  to  this  problem  for  ellipsoids  of  revolu- 
tion. Recently  Cruz-Orive  [44]  has  shown  that  if  there 
are  oblate  or  prolate  ellipsoids  with  a shape  and  size 
distribution  then  the  bivariant  three  dimensional  P.S.D. 
is  related  to  the  corresponding  bivariant  section  size 
distribution  through  a double  Abelian  integral  equation. 
Cruz-Orive  [44]  has  also  proved  that  if  there  is  a mixture 
of  oblate  and  prolate  ellipsoids  then  the  problem  does  not 
have  a unique  solution. 
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2.3.2  Distribution  of  Particle  Sizes  from  the 
Distribution  of  Section  Areas 

Johnson  [47]  and  Saltykov  [22]  have  given  algorithms 
for  the  evaluation  of  P.S.D.  from  the  distribution  of 
section  areas.  Both  Johnson  [47]  and  Saltykov  [22]  assume 
that  the  particles  are  of  spherical  shape.  Recently, 
Lewis,  Walters  and  Johnson  [48]  have  extended  and  modified 
Saltykov's  method.  Lewis,  Walters  and  Johnson  [48]  have 
presented  an  algorithm  for  a numerical  solution  to  the 
general  problem  where  the  size  class  interval  need  not  be 
constant  and  the  particles  need  to  be  of  the  same  shape 
only  within  a given  size  class.  This  is  an  excellent 
method  for  the  evaluation  of  P.S.D. , where  the  particle 
shape  may  be  a direct  function  of  its  size. 

2.3.3  Evaluation  of  P.S.D.  from  the  Distribution  of 
Linear  Intercepts 

It  is  possible  to  evaluate  the  P.S.D.  of  spheres  from 
the  distribution  of  their  linear  intercepts  by  using  the 
methods  developed  by  Spektor  [49],  Lord  and  Wills  [17]  and 
Cahn  and  Pullman  [50].  DeHoff  and  Bousquet  [51]  have  ex- 
tended Cahn  and  Pullman's  analysis  [50]  to  triaxial 
ellipsoids  of  constant  shape.  Their  results  are  sum- 
marized in  the  following  equation: 

j rij  ( £) 


(2.10) 


14 


where  rij^(£)  is  the  linear  intercept  distribution  function 
and  ny(£)  is  the  three  dimensional  P.S.D.  function.  The 
shape  factor,  K,  can  be  evaluated  by  the  method  given  by 
DeHoff  and  Bousquet  [51].  For  spherical  particles,  K is 
equal  to  2/tt.  DeHoff  and  Bousquet  [51]  have  developed  an 
expression  similar  to  this  equation  for  oriented  struc- 
tures . 

Exner  and  Lucas  [52]  have  given  a general  formula  for 
the  calculation  of  the  distribution  of  linear  intercepts 
from  a given  P.S.D.  function  for  particles  of  any  shape, 
as  long  as  all  the  particles  are  of  the  same  shape  (there 
is  a typographical  error  in  the  formula  given  by  Exner  and 
Lucas).  The  input  information  required  is  the  P.S.D. 
function  and  the  distribution  of  intercept  lengths  for  one 
particle.  The  expressions  for  the  distribution  of  inter- 
cept lengths  of  one  particle  are  available  for  the 
cube  [53],  sphere  [54],  cylinder  [55]  and  ellipsoid  [51]. 
Thus,  if  the  formula  given  by  Exner  and  Lucas  [52]  can  be 
inverted  then  it  will  give  a quite  general  solution  for 
the  evaluation  of  P.S.D.  from  the  distribution  of  linear 
intercepts . 

2.4  Evolution  of  Particle  Size  Distributions: 

The  evolution  of  microstructural  state  of 


» 


a given 
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particle  is  described  by  its  growth  path.  In  a simple 
case  where  the  geometry  of  the  particle  can  be  completely 
described  by  one  characteristic  dimension,  the  growth  path 
is  simply  the  graph  of  particle  size  against  time.  If  the 
growth  paths  of  different  particles  do  not  cross,  then  the 
number  of  particles  which  are  larger  than  a given  particle 
will  not  change  with  time.  Thus  the  number  of  particles 
larger  than  a given  particle  is  equal  to  the  particles 
that  existed  before  the  nucleation  of  this  particle.  The 
factor  which  identifies  particles  belonging  to  the  same 
growth  path,  in  a sequence  of  size  distribution  of  an 
evolving  structure,  is  that  the  N^^(R,t)  function  [see 
equation  (2.9a)]  is  constant.  This  basic  principle  was 
first  used  by  Spektor  [66]  to  evaluate  the  growth  rate  of 
carbides  during  the  coarsening  process.  Liftshitz  and 
Slyozov  [62]  and  Markworth  [63]  used  this  principle  in 
the  theoretical  analysis  of  the  coarsening  process.  The 
systematic  development  of  the  formalism  was  first  carried 
out  by  Woodhead  [57]  and  later  by  DeHoff  [58].  Woodhead ' s 
analysis  [57]  was  used  by  Mukherjee  et  al . [61]  in  their 

experimental  studies  of  particle  coarsening.*  DeHoff  [58] 

■=^Heckel  and  Degregrio  [59]  have  used  similar  analysis  to 
evaluate  the  growth  rates  of  carbide  particles  during 
spherodization  process. 
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showed  that 

n^(R,t)  = -N(t)  • (|I)^  (2.11) 

where  x is  the  nucleation  of  the  particle  whose  size  is 
R at  time  t and  N(x)  is  the  nucleation  rate  at  time  x. 

If  the  growth  rate  can  be  described  by  the  functional 

form 

(Growth  rate)  = (||)^  = f(R)  • g(t)  (2.12) 

then  equation  (2.11)  becomes 

- f(Ry.g\x)  (2-13) 

Equation  (2.13)  was  first  derived  by  DeHoff  [58]. 

Recently  Yost  [64]  has  shown  that  equation  (2.11)  is  a 
general  solution  to  the  continuity  equation  [62]  for  non- 
crossing growth  paths. 

The  global  properties  of  the  structure  can  be  evalu- 
ated from  a model  for  nucleation  rate  and  growth  rate  by 
using  equation  (2.11).  For  example,  the  volume  fraction 
of  the  transformed  phase  is  given  by  the  following  ex- 
pression : 

R 

m 

(t)  = / V(R)  n (R.t)dR  (2.14) 

ex  o 

where  V(R)  is  the  volume  of  a particle  of  size  R and  R is 

m 

the  maximum  particle  size.  For  a given  model  for  the 
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nucleation  rate  and  growth  rate,  one  can  evaluate  n^(R,t) 
at  any  time  t during  the  transformation  by  using  equa- 
tion (2.11)  and  then  equation  (2.14)  can  be  used  to  evalu- 
ate the  volume  fraction.  The  volume  fraction  evaluated  in 
this  manner  is  the  extended  volume  fraction  [the  subscript 
"ex"  on  Vy  in  equation  (2.14)  signifies  this  fact]  and  it 
can  not  be  compared  directly  with  the  experimental  values. 
The  calculation  completely  ignores  the  geometrical  im- 
pingement between  the  growing  particles.  The  impingement 
correction  becomes  increasingly  important  as  the  volume 
fraction  increases.  If  the  spatial  distribution  of  the 
particles  is  random,  then  the  following  expressions  can 
be  used  to  convert  the  extended  properties  to  their  real 
values : 


Vy  = 1 - Exp(-V^^  ) 


V 


ex 


(2.15) 

(2.16) 


ex 


M. 


1 ^^S^ex  „2 


iy=(l-Vy)[My  Sy  ] (2.16a) 

ex  (v^)^^  ex 


'V  ’ extended  volume  frac- 

ex  ex  ex 


where  V,,  . 

tion,  surface  area  and  total  curvature  calculated  from  a 

model  and  Vy,  Sy  and  My  are  their  corresponding  real 

values.  The  quantities  (v^)  and  (v„) 

S ex  ^ H ex 


will  be  defined 


18 


in  the  next  section.  Equation  (2.15)  was  first  proposed 
by  Kolmogorov  [69]  and  later  by  Avrami  [70]  and  Johnson 
and  Mehl  [71].  DeHoff  [67]  has  derived  equations  (2.16) 
and  (2.16a).  Recently  Serra  [72]  has  derived  these  equa- 
tions by  using  a Boolean  scheme.  An  alternative  deriva- 
tion of  equations  (2.15)  and  (2.16)  will  be  given  in  the 
next  chapter. 


2 . 5 Derived  Properties : 

Spektor  [65]  in  1949  showed  that  the  surface  area 
averaged  interface  velocity,  Vg,  is  given  by  the  following 
equation: 


where  v is 
ment  dS . 


//vdS 

//dS 


Sy  dt 


(2.17) 


the  local  interface  velocity  of  a surface  ele- 
Cahn  and  Hagel  [73]  rederived  this  equation  in 


1963. 


DeHoff  [67]  has  shown  that  if  the  local  interface 
velocity  is  averaged  over  the  total  curvature  of  the  sur- 
face elements,  then  the  resultant  total  curvature  averaged 
interface  velocity  is  given  by 


- _ //vHdS  _ 1 

H //HdS 


(2.18) 


In  general,  Vg  and  are  expected  to  have  different 
values . 
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In  addition  to  the  interface  velocity  averages , 
possible  to  evaluate  some  other  "average"  properties 
the  global  properties  of  a microstructure.  The  mean 
curvature,  H,  of  a phase  is  given  by  [34,35] 

5-  _ //HdS  _ 

//dS 


it  is 
from 


(2.19) 


where  H is  the  local  mean  curvature  of  a surface  element 
dS. 


For  a system  of  discrete  particles,  the  mean  inter- 
cept A of  the  particulate  phase  is  given  by  the  following 
expression  [74] : 

T 

^ = -q — (2.20) 

where  and  are  the  volume  fraction  and  surface  area 
of  the  particles. 

This  completes  the  description  of  important  stereo- 
logical  methods.  The  next  chapter  deals  with  the  theoret- 
ical analysis  developed  during  the  present  work. 


CHAPTER  III 

THEORETICAL  DEVELOPMENT 

The  geometrical  properties  of  an  evolving  micro- 
structure are  determined  by  the  following  variables : 

(1)  Starting  microstructure 

(2)  Spatial  distribution  of  nuclei 

(3)  Nucleation  rate 

(4)  Growth  rate  or  shrinkage  rate 

(5)  Particle  shape 

(6)  Growth  or  shrinkage  geometry 

An  attempt  will  now  be  made  to  answer  the  following 
questions : 

(•4)  If  a quantitative  description  of  the  above  variables 
is  given,  how  to  evaluate  the  global  properties  of  a 
microstructure? 

(B)  Given  some  global  geometrical  properties  of  micro- 
structure at  different  times  during  a phase  transforma- 
tion, how  to  evaluate  the  local  growth  rate  or  shrinkage 
rate  and  the  nucleation  rate? 

fhs  first  section  deals  with  the  analysis  of  particu- 
late structures.  More  complicated  microstructures  and  the 
impingement  problem,  etc.  will  be  analyzed  in  the  subse- 
quent sections. 
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3 • 1 Evolution  of  Particle  Size  Distributions: 
Consider  an  evolving  distribution  of  discrete  parti- 
cles in  a matrix.  The  particles  can  appear  or  disappear 
only  at  zero  size.  The  mathematical  statement  of  this 
condition  is  called  the  "continuity  equation."  Any  size 
distribution  function  must  obey  the  continuity  equation. 
The  continuity  equation,  as  it  is  given  in  the  litera- 
ture [62,64,75]  is  valid  only  if  the  growth  paths  do  not 
cross  and  it  is  given  as  follows: 


8nv(R,t) 


at 


(3.1) 


A generalized  continuity  equation  will  now  be  derived. 
Let  us  define  [see  equation  (2.9a)] 


m 

N^>(R,t)  = / n^(R,t)dR  (3.2) 

R 

where  R^  is  the  maximum  particle  size  at  time  t.  It  fol- 
lows that 


3N^^(R,t) 

( 9R )^  = -n^(R,t) 


(3.3) 


Since  the  discrete  particles  can  appear  or  disappear  only 
at  zero  size,  the  function  N^^(R,t)  is  continuous  in  both 
R and  t for  t > 0 . Thus , 


( 


3N^^(R,t) 


3R 


^R 


9N^>(R.t) 


3R 


) 


(3R) 

^3t^N 


V> 


(3.4) 
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Substituting  equation  (3.3)  into  (3.4)  and  operating  on 
the  resulting  equation  by  (-^)j-  leads  to  the  following 
equation : 


3 g 

3t  ^ 


(9R) 

^9t^N^>} 


(3.5) 


Reversing  the  order  of  differentiation  on  the  left-hand 
side  of  this  equation  and  again  substituting  equa- 
tion (3.3)  into  the  resulting  equation  gives  the  fol- 
lowing expression: 


3t  3R^^V^^’^^  ' ^yt^N 


} = 0 


V> 


(3.6) 


DeHoff  [58]  has  shown  that,  for  the  non-crossing  growth 

Thus,  for  non-crossing  growth 


paths,  (|f)^  = (II) 


3t^N 


V> 


paths,  equation  (3.6)  reduces  to  the  equation  (3.1)  How- 
ever, if  the  growth  paths  do  cross,  then  equation  (3.6) 
must  be  used.  In  order  to  solve  the  problem  of  crossing 


growth  paths,  it  is  necessary  to  formulate  (3R/3t) 


N 


in 


'V> 


terms  of  the  local  growth  rates  and  then  solve  equa- 
tion (3.6)  for  n^(R,t). 

For  a system  of  discrete  particles  in  a matrix,  one 
can  use  the  method  developed  by  Woodhead  [57]  and 


DeHoff  [58]  to  evaluate  the  growth  rate  or  shrinkage  rate, 
of  the  particles  if  the  size  distribution  is  measured  at 
different  times  during  the  evolution  of  such  a structure. 
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Fu.1T  th  eirmo  IT  6 , ori6  can  use  equation  (2.11)  to  calculate  the 
particle  size  distributions  from  a model.  These  calcu- 
luted  size  distributions  can  then  be  directly  comoared 
with  the  experimental  distributions  to  check  the  model. 
However,  equation  (2.11)  can  not  be  used  for  a dissolution 
Process.  A formalism  will  now  be  developed  for  the  pre- 
diction of  the  particle  size  distributions  during  a dis- 
solution process. 

Let  (py(^)  be  the  initial  size  distribution  function. 
Thus,  the  ^ values  represent  the  initial  particle  sizes 
which  exist  in  the  system.  Let  n^(R,t)  be  the  size  dis- 
tribution function  at  time  t during  the  dissolution  pro- 
cess. If  the  dissolution  paths  do  not  cross,  then 

R S 

m ‘^m 

/ n^(R,t)dR  = / c()^(OdC  (3.7) 

R ?[R.t]  ^ 

where  is  the  maximum  particle  size  in  the  starting 
structure.  Operating  on  the  above  equation  by  (-^)  gives 

n^(R,t)  = (f.^[C(R,t)]  (3.8) 

where  ^(R,t)  is  the  size  of  the  particle  in  the  initial 
structure,  whose  size  is  R at  time  t.  If  the  shrinkage 
rate  can  be  modeled  by  the  following  functional  form, 

(f^)^  = f(R)  • g(t)  (3.9) 

then  it  can  be  shown  that 
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.9g(R,t).  _ f(0 

9R 


IW 


(3.10) 


Combining  equations  (3.8)  and  (3.10)  gives 

<P^[^(R.t)]  • f(0 

^y(R,t)  - (3.11) 

It  is  possible  to  derive  an  expression  for  the 
annihilation  rate.  Suppose  the  particle  which  just  dis- 
appeared at  time  t had  a size  ^^(t)  in  the  initial  dis- 
tribution. For  non-crossing  dissolution  paths,  one  may 
write 


m oo  . 

/ <})y(Oci?  = -/  A(t)dt 

Co(t)  ''  t 

where  A(t)  is  the  annihilation  rate  at  time  t 
on  this  equation  by  ^ yields 

d C 


A(t)  = - 


o 


dt 


(3.12) 
Operating 

(3.13) 


Since  and  — are  completely  determined  by  the  shrink- 
age rate,  equation  (3.13)  demonstrates  that  the  annihila- 
tion rate  is  not  an  independent  variable.  One  can  write 
the  following  expression  for  the  volume  fraction  of  the 
particles : 


’m 

V^(t)  = / V(^,t)  (D^(OdC 

c„(t) 


(3.14) 


where  V(^,t)  is  the  volume  of  particle  at  time  t,  whose 
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size  was  ^ in  the  initial  distribution.  Suppose  that  the 
particles  are  of  spherical  shape  and  the  shrinkage  rate 
can  be  modeled  by  the  following  functional  form, 


^dtJ^ 


n-R' 


,n-l 


(3.15) 


where  is  the  constant  determined  by  the  kinetic  para- 
meters such  as  the  diffusion  coefficient  and  the  exponent 
n is  determined  by  the  rate  controlling  mechanism.  In- 
tegration of  equation  (3.15)  yields 

R = [c"  - (3.16) 


and 


For  spherical  particles,  combining  equations  (3.14) 
(3.16)  yields 


(3.17) 

and 


/ ^m  _ , 

V^(t)  = ^ / 4)v(OdC  (3.18) 

^o 

If  't>y(C)  is  known  then  the  parameters  n and  can  be 
evaluated  as  follow. 

(1)  Assume  a value  of  n and  calculate  V^(t)  for  different 
values  of  by  using  equation  (3.18)  Let  V^,  V^, 
etc.  be  the  experimentally  measured  values  of  the  volume 
fraction  at  times  t^^,  ^2’  • • • > stc.  One  can  calculate 
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’ ^o’  • • • > ^tc.  at  times  t^,  t2>  ■ . ■ , etc.  by  comparing 
the  experimental  values  of  volume  fraction  with  the  cal- 
culated values.  Now  plot  vs . t for  the  values  of  E 

o 

calculated  in  this  manner.  If  the  plot  is  linear,  then 

the  assumed  value  of  n is  correct  and  the  slope  of  the 

straight  line  gives  the  parameter  G . If  the  plot  of 

o o 

vs.  t is  nonlinear  then  assume  a different  value  of  n and 
repeat  the  calculations . 

Similar  procedures  can  be  developed  for  more  compli- 
cated shrinkage  rate  models.  This  algorithm  is  superior 
to  the  growth  path  analysis  because  it  requires  measure- 
ment of  only  one  particle  size  distribution.  However, 
the  technique  is  limited  to  the  dissolution  process. 

This  work  together  with  Woodhead's  [57]  and 
DeHoff's  [58]  analysis  presents  answers  to  the  questions 
(A)  and  (B)  posed  in  the  beginning  of  this  chapter  for  a 
system  of  discrete  particles,  provided  that  the  growth 
paths  do  not  cross  and  the  particle  size  distributions  can 
be  measured. 

3 . 2 Calculation  of  the  Extended  Global  Properties 
from  the  Kinetic  Models: 

Let  V(t,x),  S(t,T)  and  M(t,i)  be  the  volume,  surface 
area  and  total  curvature  of  a particle  at  time  t,  whose 
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nucleation  time  is  x.  Let  N(x)  be  the  nucleation  rate  at 
time  X.  One  can  write  the  following  expressions  for  ex- 
tended volume  fraction,  , surface  area,  S„  , and 

ex  ex 

total  curvature,  , of  a product  phase  at  time  t during 

ex 

a phase  transformation: 


t 

Vy  (t)  = / V(t,T)N(x)dx  (3.19) 

ex  o 


t 

Sy  (t)  = / S(t,x)N(x)dT  (3.20) 

ex  o 

t 

My  (t)  = / M(t,T)N(x)dx  (3.21) 

ex  o 

It  will  be  assumed  that  the  kinetic  state  of  a particle 
can  be  completely  described  by  two  size  parameters , R and 
a.  One  may  write  the  expressions  for  volume,  V,  surface 
area,  S,  and  total  curvature,  M,  of  a particle  in  terms  of 
R and  H,  if  the  quantitative  particle  shape  is  known. 
However,  almost  all  the  simple  shapes  (sphere,  cube,  cy- 
ellipsoid  of  revolution,  etc.)  and  growth  geo- 
metries can  be  treated  as  the  special  cases  of  the  fol- 
lowing model  for  which  the  volume,  surface  area  and  total 
curvature  are  given  by 


o 


£ iQ 


,2  r R nP 


If"] 

o 

R -,a 


(3.22) 


£ -.Q 


(3.23) 
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where  K^,  Kg^,  Kg^, 


(3.24) 

Kj^  3.re  the  shape  factors , 


while  and  are  the  constants  which  have  dimensions 
of  length.  The  exponents  P,  Q and  a are  constants. 

The  following  kinetics  model  will  be  assumed  for  the 
nucleation  and  growth. 

(1)  The  following  expression  gives  the  total  number  of 
particles  per  unit  volume,  Ny(t) , at  any  given  time,  t: 


N ^a+1 

yvw  - i.y  r t 


N„(t)  = N + 
o 


(3.25) 


where  , N and  a are  constants. 

'^o 

(2)  The  following  kinetic  model  gives  the  time  variation 
of  R and  £. 


(3.26) 

[Gh(t“  - 

tW)]1/S 

(3.27) 

where  and  are  the  kinetic  constants  determined  by 
the  diffusion  coefficient,  equilibrium  concentrations, 
etc.  and  the  exponents  m,  n,  co  and  6 are  constants.  Dif- 
ferent values  of  these  constants  give  different  physical 
models.  For  example,  values  of  m = 1 and  n = 2 describe 
semi-infinite  diffusion  controlled  growth. 

Combining  equations  (3.22),  (3.26)  and  (3.27)  gives 
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V(t,T) 

= (3.28) 

The  extended  volume  fraction  consists  of  two  contribu- 
tions; one  from  the  particles,  , present  at  t = 0 and 

'^o 

the  second  from  the  particles  which  nucleate  as  the  trans- 
formation proceeds  [see  equation  (3.25)] . Let  these  two 


contributions  be 

(1) 

Vy  , is  given  by 


and  V. 


(2) 

V 

ex 


ex 


V 


= N 


V 


The  first  contribution, 


V(t,o) 


ex  ■ o 

where  V(t,o)  is  the  volume  of  particles  at  time  t,  which 
were  present  at  t = 0.  Substituting  t = 0 in  equa- 
tion (3.28)  gives  V(t,o).  Thus, 


= N r r 
^PQ  ^ 

ex  o ^ 


^ n e J 


(3.29) 


where 


r = V o^~Q  /-'P/n  ,-,0/3 

^PQ  ^o  ^o  S • 


(3.29a) 


Pif f ^^Gntiating  equation  (3.25)  gives  the  nucleation  rate 
Thus,  the  nucleation  rate,  N(x),  at  time,  x,  is  given  by 


a 


N(x)  = N • x^  (3.30) 

Combining  equations  (3.19),  (3.28)  and  (3.30)  gives 

I Qco  , 111 

L"^  + o + a + 1] 


V<2)  . N Cp  D t 
ex  X » 


+ a + 1] 


(3.31) 


30 


where 


= / [1-u'^]^/’^  [l-u^^]^/^  . u""  du 
o 

(3.31a) 

Note  that  is  just  a number. 

The  total  extended  volume  fraction,  , is 

by 

given 

V ''v  '^v 

ex  ex  ex 

(3.32) 

Combining  equations  (3.29),  (3.31)  and  (3.32)  gives 


(?m  Qu 

+ n d t“-^^  . t ^ 

ex  ^ o 

(3.33) 

From  equation  (3.28),  the  volume  of  the  particle  of  maxi- 


mum  size,  V(t,o),  is  given  by 

+ 02^') 

V(t,o)  = CpQ  • t ^ 

(3.34) 

Substituting  equation  (3.34)  into  the  equation  (3.33) 
yields 


= V(t,o)  {N  + N D 

ex  '' 

(3.35) 

Similar  calculations  for  the  extended  surface  area  give 
the  following  result: 
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V ' A t”'’/"  . [N  + N D t“+l] 
ex  ^ ''o 


+ \q  t 


^ • [N^  + N Dg 


(3.36) 


where 


^P  = Gr''" 

«IQ  = Ks,  ‘o"'^  Gr'" 


(3.36a) 

(3.36b) 


D = / d-u"')^/’^  . u“  du 
^1  o 


(3.36c) 


^So  " ^ (1-u™)^/’^  (l-u“)^/^  . u“  du  (3.36d) 
2 o 

Note  that  Ap,  B^q,  Dg^  and  Dg  are  constants.  For  most  of 
the  physically  meaningful  growth  models,  either 

(1)  one  of  the  two  terms  on  the  right-hand  side  of 
equation  (3.36)  is  negligible,  or 

(2)  both  the  terms  on  the  right-hand  side  of  equa- 
tion (3.36)  have  comparable  values  but 

^1 

D . 

^2 

Thus  one  can  approximate 

/•am  Qo). 

^^ex  " ^ ].[N^  +NDgt°‘‘^^]  (3.37) 

Note  that 


/am  , Qgjv 

S(t,o)  - Apt"'’/"  + t " 


(3.38) 
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where  S(t,o)  is  the  surface  area  of  a particle  nucleated 
at  t = 0.  Combining  equation  (3.37)  and  (3.38)  gives 

= S(t,o)  [N„  + N D (3.39) 

ex  ^ 

where  Dg  is  equal  to  Dg  if  the  second  term  in  equa- 
tion (3.36)  is  negligible  and  it  is  equal  to  D„  if  the 
first  term  can  be  neglected.  Dg  is  approximately  equal  to 
either  Dg^  or  Dg  if  both  the  terms  on  the  right-hand 
side  of  equation  (3.36)  have  comparable  values. 

Calculations  for  the  extended  total  curvature  yield 
the  following  result: 

am 

“v  tNy  + ND 

ex  '^o  “l 

Qo) 

+ E t^  [N  +ND  t“+^]  (3.40) 

^ o 2 

where 


= K,,  £ 


1-Q  pQ/6 


"Q  "M2  ^o  "L 


= / [1-u'^]^/''  u""  du 

1 o 

= / [l-u"^]Q/e  u“  du 

2 o 


(3.40a) 

(3.40b) 

(3.40c) 


(3.40d) 
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Note  that  E^,  Eq,  and  are  constants.  Again  using 

the  similar  arguments  as  those  used  in  the  case  of  ex- 
tended surface  area,  one  can  approximate 

= M(t,0)  [N  + N D t""+^]  (3.41) 

ex  ''o 

where  M(t,0)  is  the  total  curvature  of  a particle  nucle- 
ated at  t = 0 . 

The  various  special  cases  of  this  model  will  be  ao- 
plied  to  real  systems  in  the  next  three  chapters.  It  is 
possible  to  develop  more  sophisticated  models  by  assuming 
some  other  functional  forms  for  the  nucleation  rate  and 
growth  rate.  For  example,  if  the  total  number  of  parti- 
cles per  unit  volume,  N^(t) , has  the  following  functional 
form. 


j.  N ,_a+l  , 

INy^  + t + 

then  the  expression  for  the  extended 
given  by 


C .Y+1 
Ty+IT 

volume  fraction 


(3.42) 

is 


Vy  = V(t,0)  [Ny  + N D t""^^  + C D„  t^^^]  (3.43) 

ex  o ^ 

where 


= / (1-u™)^'^’^  (l-u“)^/^  u^  du  (3.43a 

o 

and  C and  y are  constants.  However,  the  present  model  is 
flexible  enough  to  be  valid  for  most  of  the  cases. 
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3 . 3 Evaluation  of  the  Nucleation  Rate  and 
Growth  Rate  from  the  Global  Propertied 

1^  this  section,  an  attempt  will  be  made  to  answer 
the  question  (B)  posed  in  the  beginning  of  this  chapter. 
It  may  be  said  that  in  most  of  the  physically  meaningful 
cases  a given  microstructure  at  time  t during  a phase 
transformation  contains  a complete  history  of  the  fol- 
lowing variables: 

(1)  Starting  microstructure 

(2)  Nucleation  rate 

(3)  Growth  rate 

(4)  Spatial  distribution  of  nuclei 

(5)  Particle  shape  and  Growth  geometry 

this  is  so,  then,  in  principle,  it  should  be  possible 
to  invert  this  problem.  Thus  it  should  be  possible  to 
evaluate  the  nucleation  and  growth  rate,  etc.  if  all  the 
global  properties  of  the  growing  phase  are  measured  at 
different  times  as  the  transformation  proceeds.  However, 
at  present  the  stereological  techniques  have  not  been  de- 
veloped to  measure  al 1 the  global  properties  of  the 
growing  phase  from  the  plane  of  polish.  Thus  some 
simplifying  assumptions  regarding  these  five  variables 
may  be  needed  to  solve  the  problem. 

A technique  will  now  be  developed  to  evaluate  the 
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nucleation  rate  and  growth  rate  from  the  measured  values 
of  the  volume  fraction,  surface  area  and  total  curvature. 
In  order  to  make  the  problem  mathematically  tractable  the 
following  assumptions  will  be  made. 

(1)  It  is  possible  to  convert  the  measured  values  of  the 
global  properties  to  their  corresponding  extended  values. 

(2)  The  growth  rate  of  a given  particle  is  independent 

of  the  nucleation  and  growth  of  the  other  particles.  This 
implies  that  the  growth  rate  of  a given  particle  depends 
only  on  the  amount  of  time  for  which  the  particle  is 
present  in  the  system.  If  a particle  nucleated  at  time, 

T,  then  its  growth  rate,  at  any  time  t during  the  trans- 
formation, depends  only  on  (t-i) . It  follows  that,  the 
volume,  surface  area  and  total  curvature  of  the  particle 
will  be  functions  of  (t-x)  only.  Thus, 

V(t,T)  = V(t-x)  (3.44) 

S(t,x)  = S(t-x)  (3.45) 

M(t,x)  = M(t-x)  (3.46) 

Note  that  many  theoretical  models  for  growth  assume  this 
situation  [3,6,8-11]. 

Combining  equations  (3.19)  and  (3.44)  gives 

t 

Vy  (t)  = / V(t-x)N(x)dx  (3.47) 

ex  o 


Similarly , 
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t 

(t)  = / S(t-T)N(T)dT  (3.48) 

ex  o 

and 

t 

(t)  = / M(t-T)N(x)dT  (3.49) 

ex  o 

The  integrals  on  the  right-hand  side  of  these  equations 
are  called  "convolution  integrals."  If  the  particle  shape 
is  known  or  at  least  can  be  modeled,  then  the  functions 
_V(t-T) , S(t-T)  and  M(t-x)  are  not  independent  of  each__ 
other.  For  a given  particle  shape,  they  are  completely 
determined  by  the  growth  rate.  In  such  a case,  there 
are  two  unknown  functions  (nucleation  rate  and  growth 
rate)  in  the  above  three  integral  equations.  If  one  can 
solve  these  integral  equations,  then  both  the  nucleation 
rate  and  growth  rate  can  be  obtained. 

It  is  possible  to  solve  the  integral  equations  (3.47) 
to  (3.49)  by  using  Laplace  transform  techniques.  A 
Laplace  transform  [80]  of  a function,  ip(t),  is  defined  as 

00 

L[4^(t)]  = / e"^^  i|;(t)dt  (3.50) 

o 

The  transform,  L[t[;(t)],  is  a function  of  the  parameter, 

S,  only.  An  inverse,  L of  a Laplace  transform  [81]  is 
defined  as  follows: 


L ^ [L{i|;(t)  } ] = (t) 


(3.51) 
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The  Laplace  transform  of  a convolution  integral  has  the 
following  interesting  properties  [82]: 
t 

L{/  V(t-T>N(T)dT}  = L{V(t)}  • L{N(t)}  (3.52) 

o 

Taking  the  Laplace  transform  of  both  sides  of  equa- 
tion (3.47)  gives 


t 

L{V^  (t)}  = L{/  V(t-x)N(T)dT}  (3.53) 

ex  o 

Combining  equations  (3.52)  and  (3.53)  gives 

L{V^  (t)}  = L{V(t)}  • L{N(t)}  (3.54) 

ex 

Rearranging  the  terms  in  equation  (3.54)  and  taking  the 
inverse , 


N(t) 


L[V^  (t)] 


{ ex  . 
^ L[V(t)j  ^ 


(3.55) 


Similarly  equations  (3.48)  and  (3.49)  give  the  following 
result : 


L[S^  (t)j 
N(t)  = L~^{- 


L[S(t)] 


L[M^  (t)j 


(3.56) 

(3.57) 


Thus  if  the  growth  rate  is  known  then  the  nucleation  rate 
can  be  evaluated  by  using  any  of  these  three  equations. 

Combining  the  above  equations  for  N(t) , two  at  a time 
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and  using  the  convolution  theorem  [82] , leads  to  the  fol- 
lowing expressions 


/ V(t-T)S^  (T)di  = / S(t-x)V„  (x)d 


o 


ex 


V 


(3.58) 


ex 


/ S(t-x)M^  (x)dx  = / M(t-x)S„  (T)d 


o 


ex 


V 


ex 


t t 

/ M(t-x)V^  (x)dx  = / V(t-x)M^  (x)dx 
o ex  o ''ex 


(3.59) 


(3.60) 


A careful  examination  of  these  equations  leads  to  the  fol- 
lowing conclusions : 

(1)  Only  two  of  these  three  equations  are  independent. 

(2)  The  nucleation  rate  does  not  appear  explicitly  in 

these  equations.  The  quantities  V„  , S„  and  M,, 

V V V are 

ex  ex  ex 

obtained  from  the  experimental  data,  while  V(t-x) , S(t-x) 
and  M(t-x)  are  determined  by  the  growth  rate.  Thus  the 
quantitative  evaluation  of  the  local  growth  rate  is  pos- 
sible without  any  knowledge  of  the  nucleation  behavior. 

If  the  particle  shape  is  known,  then  V(t-x) , S(t-x) 
and  M(t-x)  are  not  independent  unknown  functions.  For 
example,  if  the  particles  have  a spherical  shape  then 


S(t-T)  = 2^[M(t-x)]2 


(3.61) 


and 


V(t-x)  = — L^[M(t-x)]^  (3.62) 

48tt^ 
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In  such  a case,  substituting  equations  (3.61)  and  (3.62) 
in  equations  (3.58)  and  (3.59)  gives 

^3  ^2 

(T)di  = 12TT/[M(t-T)  ]^V  (x)dT  (3.63) 

o ex  o ''ex 

t „ t 

/[M(t-x)]  My  (x)dx  = 4fT/M(t-x)Sy  (x)dx  (3.64) 
o ex  o ''ex 

The  third  equation  given  by  equation  (3.60)  is  not  an 
independent  equation.  There  is  only  one  unknown  function, 
M(t-x)  , in  the  above  two  equations.  One  can  now  solve 
equation  (3.63)  by  numerical  techniques  [84]  to  get  the 
values  of  M(t-x)  for  the  different  values  of  the  argument. 
The  solution  can  be  checked  by  using  equation  (3.64).  Al- 
ternatively, if  several  models  are  available  for  the 
growth  process,  then  they  can  be  checked  by  using  these 
equations,  i^ote  that  this  technique  does  not  require  any 
icnowledge  or  assumption  regarding  the  nucleation  behavior. 
Furthermore,  the  data  processing  involved  in  the  equa- 
tions (3.63)  and  (3.64)  is  the  integration  of  the  experi- 
mental data.  Thus,  the  processing  errors  are  expected  to 
be  small.  This  approach  will  be  used  in  the  fourth 
chapter  to  deduce  the  growth  rate  and  the  nucleation  rate 


of  austenite  during  an  inverse  eutectoid  transformation. 
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3 . 4 Growth  Geometry: 

Various  growth  geometries  will  be  described  in  this 
section.  The  importance  of  the  growth  geometry  in  the 
interpretation  of  experimental  data  will  be  demonstrated. 
In  what  follows,  the  growth  geometries  will  be  described 
for  a cylindrical  particle  of  radius  R and  length  £.  The 
description  can  be  extended  to  some  other  particle  shapes . 

3.4.1  Shape  Preserving  Growth 

During  shape  preserving  growth,  the  quantitative 
particle  shape  does  not  change  with  time.  It  follows  that 
spherical  growth  is  a special  case  of  shape  preserving 
growth.  If  the  growth  of  a cylindrical  particle  occurs 
in  such  a way  that  (^i/Ry  does_not  change  with  time.,  then 
it  can  be  called  a shape  preserving  growth.  The  model 
developed  in  section  3.2  will  now  be  analyzed  for  the 
shape  preserving  growth  of  cylinders. 

^ = Constant  = Ag  (3.65) 

Comparing  equations  (3.26),  (3.27)  and  (3.65)  gives 

m = w;  n = 6 and  G^  = Ag  Gj^  (3.66) 

Furthermore,  since  the  particles  have  a cylindrical  shape, 
the  shape  factors  in  equations  (3.22)  to  (3.24)  take  the 
following  values: 
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Ky  = TT 

{ 


(3.67) 


The  exponents  in  equations  (2.22)  to  (2.24)  take  the 
following  values: 


{P  = 2;  Q = a = 1}  (3.68) 

Substituting  these  values  in  the  equation  (3.33)  gives 

Vy  = Cy-,  t^"'/’^[Ny  + N Dy  (3.69) 

ex 

where  C21  = (Ag)^^^  • tt 

and  Dy  = /[1-u"^]^/^  • u“  du 
o 

Similarly  for  this  particular  case,  equations  (3.36)  and 
(3.40)  give  the  following  results: 

S^  = {2tt[1+Ac]  • + N D„  t“'^^](3.70) 

ex  a K ^o  ^ 

Note  that  Dg  = Dg  = Dg  for  the  shape  preserving  growth. 

My  = {[7T^+  Ag]G^/"'}t'”/’^[Ny  + N t“+^]  (3.71) 

One  can  now  calculate  (v^.)  and  (v,t)  by  using  equa- 

O 0X.  1*1  0X 

tions  (3.69)  to  (3.71).  The  behavior  of  (v„)  and  (v„) 

O 0X  ii  0X 

is  determined  by  all  the  parameters  of  the  model  (Gj^,  m, 
n,  etc.).  Thus  in  this  case  the  interface  velocities  may 
be  used  to  get  information  about  the  growth  rate;  it  will 
now  be  shown  that  this  is  not  always  true. 
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3.4.2  Two  Dimensional  Growth 

If  a cylindrical  particle  grows  in  such  a manner 
that  its  radius  R increases  with  time,  while  its  length 
remains  constant,  then  such  a growth  behavior  is  called 
"two  dimensional  growth."  The  concept  of  two  dimensional 
growth  can  be  extended  to  some  other  particle  shapes. 

Let  us  denote  the  constant  length  of  the  cylindrical 
particle  by  One  can  write  the  following  expressions: 

V(t,T)  = [tt£^]  [R(t,T)]^  (3.72) 

S(t,x)  = 27T[R(t,T)]^  + 2tt£^  • R(t,x)  (3.73) 

Since  the  radius  of  the  particle  increases  with  time  while 
its  length  remains  constant,  soon  after  nucleation  the 
first  term  on  the  right-hand  side  of  the  equation  (3.73) 
will  be  much  larger  than  the  second  term.  Thus, 

S(t,x)  = 2TT[R(t,x)]^  (3.74) 

Substituting  equations  (3.72)  and  (3.74)  into  the  equa- 
tions (3.19)  and  (3.20)  gives 

t 2 • 

Vy  (t)  = TT£^/[R(t,x)]^  N(x)dx  (3.75) 

ex  o 

and 

t 2 • 

S^  (t)  = 2TT/[R(t,x)  N(x)dx  (3.76) 

^ex  o 

Substituting  equation  (3.76)  into  (3.75)  gives 
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Vy  (t)  = (-^)  Sy  (t)  (3.77) 

ex  ex 

Since  does  not  change  with  time,  equation  (3.77)  pre- 
dicts a linear  relationship  between  the  extended  volume 
fraction  and  surface  area.  Equation  (3.77)  is  a direct 
consequence  of  the  assumed  two  dimensional  growth  geo- 
metry. No  assumptions  regarding  the  nucleation  rate  or 
growth  rate  were  made  in  deriving  equation  (3.77).* 

The  extended  surface  area  averaged  interface  velo- 
city, (Vg)g^,  is  given  by 


d V 


V 


ex 


ex 


V 


dt 


(3.77a) 


ex 


For  two  dimensional  growth,  substituting  equation  (3.77) 
into  equation  (3.77a)  yields 

d In 

(3.78) 


£ 


N O 

‘'^S^ex  2 


dt 


The  behavior  of  (^g)gj^  will  now  be  analyzed  for  the  model 
developed  in  section  3.2.  For  two  dimensional  growth  of 
cylindrical  particles,  the  parameters  P,  Q and  a in  the 
equations  (3.22)  to  (3.24)  take  the  following  values: 

{P  = 2;  Q = 0;  a = 1}  (3.79) 


*Note  that  for  a two  dimensional  growth,  equation  (2.58) 
becomes  redundant. 
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For  these  values  of  the  parameters,  equation  (3.33)  takes 
the  following  form: 


V 


V. 


(t)  = + N (3.80) 


’V 


V 


ex  ' ■■  ’o 

Substituting  equation  (3.80)  into  (3.78)  gives 

a 


_ £ m , 

('’S>ex  - + 


nND^(a+l)t 

^ m(N^  +ND^t“'''^) 
'^o 


J 


If  Ny  =0,  then 
'^o 


£ m 


X.  Ill  / , 1 \ n 

• (^)Cl  + ■ i 


"S^ex  2n  ' m 

If  N = 0 (i.e.  site  saturation),  then 


(3.81) 


(3.82) 


_ £ m , 

(«S>ex  = • E <3.83) 

The  coefficients  of  1/t  in  equations  (3.82)  and  (3.83) 
do  not  depend  on  temperature.  These  equations  predict 
that  X^g)^^  should  be  independent  of  temperature,-- 
Furthermore,  the  hyperbolic  time  dependence  of  (v^)  is 
not  affected  by  the  parameters  of  the  growth  mode  (i.e. 
m and  n)  . 

This  analysis  of  two  dimensional  growth  will  be  used 
in  the  next  chapter  for  the  interpretation  of  some  experi- 
mental data. 
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3.4.3  One  Dimensional  Growth 

If  a cylindrical  particle  grows  in  such  a manner  that 
its  length,  Si,  increases  with  time  while  its  radius  re- 
mains constant,  then  such  a behavior  is  called  "one 
dimensional  growth."  This  concept  can  be  extended  to  some 
other  particle  shapes.  Let  us  denote  the  constant  radius 
of  the  cylindrical  particle  by  . One  can  write  the 
following  expressions: 


V(t,x)  = [ttR^]  £(t,x) 

(3.84) 

S(t,x)  = [2ttR^]  £(t,x) 

(3.85) 

M(t,x)  = IT  • £(t,x) 

Combining  equations  (3.84)  and  (3.19)  gives 

(3.86) 

2 ^ 

(t)  = ttR^  / Jl(t,x)N(x)dx 
ex  o 

Similarly , 

(3.87) 

t 

(t)  = 2ttR  / £(t,x)N(x)dx 
ex  o 

and 

(3.88) 

t 

My  (t)  = IT  / £(t  ,x)N(x)dx 
'^ex  o 

Combining  equations  (3.87)  and  (3.89)  gives 

R . 

Vy  (t)  = (^)  S (t)  = R„  My  (t) 

ex  ex  ex 

(3.89) 

(3.90) 

(3.90) 
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Since  does  not  change  with  time,  equation  (3.90)  pre- 
dicts a linear  relationship  between  the  extended  volume 
fraction,  surface  area  and  total  curvature.  Again,  the 
derivation  of  equations  (3.90)  do  not  require  any 
assumptions  regarding  the  nucleation  rate  or  growth  rate 
Combining  equations  (3.78)  and  (3.90)  gives 


^^S^ex 


R d In  Vy 
o ex 


(3.91) 


2 dt 

The  extended  total  curvature  averaged  interface  velocity, 
^^H^ex’  given  by 


d S 


(v„) 


V 


ex 


H''  ex  2 M 


V 


dt 


ex 


Substituting  equation  (3.90)  into  (3.92)  gives 


(v„)  = R 

^ H ex 


o 


d In  S 

3T 


V 


ex 


-=  R 


o 


d In  V 

dF 


V 


ex 


(3.92) 


(3.93) 


or 

^^H^ex  = 2(Vg)ex  (3.94) 

The  behavior  of  (Vg)^^  and  (v^)^^^  will  now  be  analyzed  for 
the  model  developed  in  section  3.2.  For  one  dimensional 
growth , 

{P  = 0;  Q = 1;  a = 0}  (3.95) 

Substituting  these  values  of  P , Q and  a together  with 
shape  factors  for  the  cylindrical  [equation  (3.67)]  in 
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equation  (3.33)  gives 

Vy  (t)  = [ttR^gJ;/^]  . t*^/^  . [N 


ex 


V + N Dy  t 
o 


a+1 


](3.96) 


Substituting  equation  (3.96)  in  (3.91)  gives 

1 


_ R CO 


a 


+ 


R^ND^(a+l)t 

^ 2(N^  +ND^t“'^^) 

o 


(3.97) 


Combining  equations  (3.94)  and  (3.97)  yields 

R c 
o 

'S'ex 


ex 


- (v=).„  = <^)  • i + (3.98) 


2(N^  +ND^t 
o 


a 
a+1 


If  = 0, 
o 


ex 


R 


2 ^^S^ex  " + a + 1)] 


1 

t 


(3.99) 


If  N = 0 (i.e.  site  saturation), 
— 2 = 


R (0 


1 

t 


(3.100) 


'S'' ex  *•  23 

The  quantities  in  the  square  brackets  in  equations  (3.99) 
and  (3.100)  do  not  depend  on  the  temperature.  These  equa- 
tions predict  that  arid  should  be  independent 

of  temperature.  Furthermore,  the  hyperbolic  time  depen- 
dence of  (’^g)gj^  and  (’^jj)gj.  is  not  affected  by  the  para- 
meters of  the  growth  model  (i.e.  Gj^,  to  and  B)  . 


3 . 5 The  Impingement  Problem: 


Imagine  a structure  calculated  from  some  model  based 
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upon  an  assumed  nucleation  rate  and  growth  rate  (for 
example,  the  model  in  section  3.2).  If  the  standard 
methods  for  calculating  the  global  properties  of  the 
structure  are  used,  such  calculations  will  over  estimate 
the  corresponding  properties  in  the  real  structure.  The 
magnitude  of  this  over  estimation  increases  with  the  in- 
crease in  the  volume  fraction  of  the  growing  phase.  The 
over  estimates  result  because  the  geometrical  impingement 
of  the  growing  particles  is  not  taken  into  account.  As 
mentioned  in  the  second  chapter,  the  modeled  structure  in 
which  the  particles  may  grow  through  one  another  is  called 
the  "extended  structure."  If  the  particles  are  distri- 
buted randomly  in  the  structure,  then  equations  (2.15)  and 
(2.16)  give  relations  between  some  real  and  extended 
global  properties.  A simple  derivation  of  these  equations 
will  be  first  given. 

It  will  be  assumed  that  the  particle  centers  are 
distributed  randomly  in  the  matrix.  For  such  a case, 
Chandrasekhar  [76]  and  Hertz  [77]  have  given  the  following 
expression : 

(|)(x)  = 4ttx^  Exp{-^  (3.101) 

where  4>(x)dx  is  the  probability  that  the  first  nearest 
neighbor  of  a given  particle  center  is  in  the  spherical 
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shell  of  radii  x and  (x  + dx)  and  is  the  total  number 

of  particles  per  unit  volume. 

The  derivation  will  be  first  given  for  a simple  case 

where  all  the  particles  are  spheres  of  radius  R(t) . Let 

dVy  and  dV„  be  the  changes  in  the  real  and  extended 
'^ex 

volume  fractions  during  a small  increment  in  time,  dt. 

Focus  upon  a small  part,  6(dV^  ),  of  an  incremental 

2 ^ex 

layer,  (4ttR  dR)  , on  one  particle.  Now,  6(dVy  ) will  con- 

''ex 

tribute  to  the  real  volume  only  if  there  is  no  particle 

center  in  a sphere  of  radius  R(t)  around  it.  The  pro- 

R(t) 

bability  of  this  event  is  [1  - / (})(x)dx].  This  pro- 

o 

bability  is  the  same  for  all  such  volume  elements  in  the 
structure  at  a given  time  t and  it  gives  fraction  of  such 
elements  which  contribute  to  the  real  structure.  Hence, 

R(t) 

Z 6(dV^)  = [1  - / c!)(x)dx]  Z ) 

o ''ex 


or 


R(t) 

dVy  = [1  - / c()(x)dx]  • dV 


V 


(3.102) 


ex 


Substituting  equation  (3.101)  into  (3.102)  and  invoking 
the  fact  that 
lowing  result: 


the  fact  that  V„  is  equal  to  ^[R(t)]^N^  gives  the  fol- 

^ex  ^ ^ 


-V 


V 


dVv  = 


ex 


dV. 


V 


ex 


e 


(3.103) 
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Integration  of  this  equation  gives 

= 1 - Exp(-V^  ) 
ex 

This  equation  is  identical  to  equation  (2.15)  proposed  by 

Kolmogorov  [69]  and  Avrami  et  al . [70] . 

Consider  a small  elemental  area,  6 (S  ),  on  the  ex- 

ex  ’ 

tended  surface  of  any  particle.  This  element  will  con- 
tribute to  the  real  surface  area  only  if  there  is  no 

particle  center  in  a sphere  of  radius  R(t)  around  it.  The 

R(t) 

probability  of  this  event  is  again  [1  - / (})(x)dx].  This 

o 

probability  is  equal  to  the  fraction  of  such  elements 
which  contribute  to  the  real  surface  area.  Thus 


R(t) 

I 6(S)  = [1  - / ([)(x)dx]  Z 6(S) 

GX 


or 


R(t) 

Sy  = [1  - / 4>(x)dx] 

o ex 


But , 


R(t)  ~^V 

[1  - / (J)(x)dx]  = e = (1-V„) 

o 


Hence , 

Sy  = Sy 

ex 


This  equation  is  identical  to  equation  (2.16). 


51 


Now  let  us  consider  the  case  where  all  the  particles 
are  not  of  same  size  in  the  extended  structure.  Let 
E(R,t)  be  the  extended  particle  size  distribution  function 
of  the  spherical  particles.  Let  us  divide  this  distribu- 
tion function  into  equally  spaced  histograms  and  let 
AR  be  the  size  class  interval.  Note  that 

^e 

V = £im  { Z ^ R?  E(R. , t)AR}  (3.104) 

''ex  N i=l  ^ ^ 

e 

AR->-o 

Again,  consider  an  element,  6(dV^  ),  on  any  particle. 

''ex 

The  probability,  , that  this  volume  element  will  contri- 
bute to  the  real  structure  is  a product  of  N^  pro- 
babilities. The  probability,  P , that  there  is  no  parti- 
cle  of  radius  R^  inside  a sphere  of  radius  R^  around  this 
volume  element  is  given  by 

Ri 

P^  = [1  - / 4)j^(x)dx] 

i o 

where 


E(R^,t)AR 


Exp{ 


4itx' 

■~T“ 


E(R^,t)AR} 


Thus 


4uR? 

E(R^,t)AR 


52 


The  probability,  , that  the  incremental  voliame  element, 

6(dVy  ),  will  contribute  to  the  real  structure  is  given 
ex 
by 

e e 4ttR. 

Pr  = P^_  = Exp{ E(R^,t)AR} 


or 


N , 
e 4ttR  . 


P = Exp{-  E — ^ E(R.,t)AR} 

^ i=l  1 


As  and  AR->o,  the  summation  is  exactly  equal  to  V. 

[see  equation  (3.104)].  Thus, 


V 


ex 


P^  = Exp(-V^  ) 
ex 


(3.105) 


Again,  the  probability,  P^ , is  equal  to  the  fraction  of 

the  incremental  elements,  6(dV^  ),  which  contribute  to 

^ ex 

the  real  structure.  Hence, 

E 6(dV^)  = {Exp(-V^  )}  E 6(dV^  ) 


ex 


ex 


or 

dV^  . e dV 

ex 

The  integration  of  this  equation  gives  equation  (2.15). 
To  derive  equation  (2.16),  consider  a small  surface  ele- 
ment, 6(S  ),  on  the  extended  surface  of  any  particle. 

The  probability  that  this  element  will  contribute  to  the 
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real  structure  is  exactly  equal  to  . Thus, 

Z 6(S)  = {Exp(-V  )}  E 6(S^) 

or 

Sy  = {Exp(-V^  )} 

ex  ex 

or 

Sy  = (1  - V^) 

ex 

The  approach  can  be  generalized  for  convex  particles  of 
other  shapes  if  their  orientation  is  also  random. 

3.5.1  Non-random  Impingement 

Let  us  assume  that  all  the  particles  are  spheres  of 
radius  R(t) , but  the  distribution  of  their  particle  cen- 
ters is  non-random.  Consider  an  incremental  element, 

5 (dV.^  ),  on  the  surface  of  a particle.  Let  the  spatial 

coordinates  of  this  element  be  (u,  v,  w) . The  probability 
that  there  is  no  particle  center  in  a sphere  of  radius  x 
and  there  is  at  least  one  particle  in  a spherical  shell  of 
radii  x and  (x  + dx)  is  given  by 

X 

(})(x)dx  = [1  - /(})(x)dx]  P (x)dx  (3.106) 

U,  V Vv 

o 

where  the  first  term  on  the  right-hand  side  of  this  equa- 
tion is  the  probability  that  there  is  no  particle  center 
in  a sphere  of  radius  x.  The  second  term, 
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the  probability  that  there  is  at  least  one  particle  center 
in  the  spherical  shell  of  radii  x and  (x  + dx) . The  sub- 
script "uvw”  signifies  that  this  probability  will  depend 
on  the  spatial  coordinates  (u,  v,  w)  of  the  incremental 

element,  6(dV„  ).  Let  V be  the  total  volume  of  the 

''ex 

sample.  An  average  value,  P(x)  , of  can  now  be 

defined  as  follows: 

P(x)  = ^ ///  (x)dudvdw  (3.107) 

''o 

Since  we  are  interested  only  in  the  average  properties, 
a "spatial  average"  value  of  <j)(x)  can  be  evaluated  by 
using  the  following  equation: 

X _ 

(|)(x)  = [1  - /4>(x)dx]  P(x)  (3.108) 

o 

The  solution  to  this  integral  equation  gives  the  following 
result : 

_ x_ 

4>(x)  = P(x)  Exp{-/P(x)dx}  (3.109) 

o 

The  probability  that  there  is  no  particle  in  a sphere  of 
radius  R(t)  is  given  by 

R(t)_ 

pr  = Exp{-  / P(x)dx}  (3.110) 

o 

Thus,  if  all  the  particles  are  spheres  of  radius  R(t) , 
then 
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2 

dV^  = 47T[R(t)j  N^[Exp{-  / P(x)  dx}]dR  (3.111) 

o 

or 

R(t)  . R(t)_ 

= / 47T[R(t)]  Exp[-  / P(x)dx]dR(t)  (3.112) 

o o 

The  development  can  be  generalized  if  there  is  a distribu- 
tion of  spheres  in  the  extended  space.  Thus,  if  quanti- 
tative information  about  P(x)  is  available,  then  the  non- 
random  impingement  problem  can  be  solved.  However, 
physically  meaningful  models  for  P(x)  must  be  developed 
for  the  further  analysis  of  this  approach. 

3.5.2  A Phenomenological  Approach  to  Non-Random 
Impingement 

A simple  approach  will  now  be  presented  which  is 
extremely  useful  for  two  dimensional  and  one  dimensional 
growth.  It  may  be  conjectured  that  the  physical  factors 
responsible  for  these  growth  geometries  may  also  give  rise 
to  non-random  distribution  of  particle  centers  in  space. 
For  two  dimensional  or  one  dimensional  growth  [see  equa- 
tions (3.77)  and  (3.90)] , 

Vv  = K • S„  (3.113) 

ex  ex 

where  K is  a geometrical  constant. 

It  will  be  assumed  that 
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dVv  = nVy)  dv^ 


(3.114) 


ex 


where  'P  (Vy)  is  the  impingement  function  to  be  evaluated, 
Now, 


V 


V, 


ex 


' o 


(3.115) 


Using  the  same  arguments  as  used  earlier, 

Sx 


"V 


ex 


'V 

¥(V^ 


(3.116) 


Combining  equations  (3.113),  (3.115)  and  (3.116)  yields 

V, 


V dV 


/ 


V 


= K 


'V 


(3.U7) 


The  solution  to  this  integral  equation  is  as  follows; 

dV, 


■/ 


V 


KS 


^(Vy)  = K Sy  e 


V 


+ B (3.118) 

where  B is  a constant  and  can  be  evaluated  as  follows: 


dV 


-/ 


V 


KS 


B = - lim  {K  Srr  e 
''V*! 


V 
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3 . 6 Conclusions : 

(1)  A generalized  continuity  equation  has  been  de- 
rived and  an  approach  to  the  solution  of  the  problem  of 
crossing  growth  paths  is  suggested. 

(2)  An  expression  is  derived  for  the  particle  size 
distribution  function, during  a dissolution  process,  in 
terms  of  the  initial  particle  size  distribution  function 
and  shrinkage  rate. 

(3)  An  algorithm  is  developed  for  the  evaluation  of 
shrinkage  rate  during  a dissolution  process . The  input 
information  required  is  the  initial  particle  size  distri- 
bution function  and  the  values  of  the  global  properties  of 
the  dissolving  phase  at  various  times  during  the  dissolu- 
tion process . 

(4)  A systematic  approach  has  been  developed  for 
the  evaluation  of  the  global  properties  of  a growing 
phase  from  a given  description  of  the  local  variables. 

A flexible  but  simple  model  is  developed  for  the  calcula- 
tion of  global  properties . 

(5)  A fairly  general  technique  has  been  developed 
to  deduce  the  nucleation  rate  and  growth  rate  from  the 
time  variation  of  the  global  properties. 

(6)  Effects  of  various  growth  geometries  on  the 
kinematics  of  phase  transformation  are  explored  in  detail. 
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(7)  A simple  derivation  of  the  impingement  equations 
is  given  for  the  case  of  random  impingement. 

(8)  A solution  to  the  general  non-random  impingement 
problem  is  suggested.  A phenomenological  approach  is 
developed  for  the  evaluation  of  the  non-random  impingement 
function  for  two  dimensional  and  one  dimensional  growth. 


CHAPTER  IV 
RECRYSTALLIZATION 

4 . I Introduction : 

The  purpose  of  this  chapter  and  the  next  two  chapters 
is  to  demonstrate  the  applicability  of  the  analysis  de- 
veloped in  the  previous  chapters.  In  this  chapter,  the 
data  of  some  other  research  workers  [79,92]  on  the  kine- 
tics of  the  recrystallization  process  will  be  reanalyzed 
and  reinterpreted  by  using  the  present  approach.  It  will 
be  shown  that  some  apparent  "paradoxes"  which  exist  in  the 
literature  can  be  resolved  by  using  some  simple  explana- 
tions based  only  on  the  geometry  of  the  microstructural 
evolution . 

General  reviews  of  the  detailed  treatments  on  the 
topic  of  recrystallization  can  be  found  in  different 
works  [95-97].  Nucleation  in  recrystallization  is  de- 
fined as  the  formation  of  a strain-free  region  whose 
boundary  can  transform  the  strained  matrix  into  a strain- 
free  material  as  it  moves . Two  mechanisms  are  identified 
for  the  process;  one  involving  the  existing  grain  bounda- 
ries (strain  induced  boundary  migration) , and  the  second 
involving  the  regions  of  sharp  lattice  curvature  (sub- 
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grain  growth) . Theoretical  models  for  the  kinetics  of 
grain  boundary  migration  have  been  developed  [98-101]. 

In  experimental  studies,  the  kinetics  of  recrystalli- 
zation may  be  monitored  by  measuring  the  global  properties 
of  the  microstructure  during  the  process.  Speich  and 
Fisher  [79]  have  studied  the  recrystallization  of  3%7o 
silicon  steel  by  using  quantitative  microscopic  techni- 
ques. The  experimental  data  of  Speich  and  Fisher  [79] 
will  now  be  reanalyzed. 

4 . 2 The  Experimental  Work  of  Speich  and  Fisher  [79] : 

The  experiments  were  performed  on  a 3.277o  silicon 
steel  with  0.0197o  carbon.  The  grain  size  of  the  starting 
material  was  about  fifty  microns.  The  material  was  cold 
rolled  607o  to  produce  a strained  matrix.  The  recyrstal- 
lization  of  this  strained  matrix  was  studied  at  various 
temperatures  between  500°C  to  1000°C.  A laser  heating 
technique  was  used  for  short  annealing  times  (0.002  sec 
to  0.1  sec  used  at  high  temperatures).  The  volume  frac- 
tion, Vy , and  the  surface  area,  Sy,  of  the  recrystallized 
regions  was  measured  by  using  the  standard  stereological 
techniques . 
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4.2.1  Experimental  Observations 

(1)  The  recyrstallized  grain  size  was  found  to  be 
independent  of  the  annealing  temperature.  Note  that  the 
experiments  were  performed  at  various  temperatures  ranging 
from  500°C  to  1000°C.  The  recrystallized  grain  size  was 
about  11  microns. 

(2)  The  relationship  between  the  volume  fraction, 

V^,  and  surface  area,  S^,  could  be  described  for  all  the 
times  and  temperatures  studied  by  a single  empirical  equa- 
tion : 


^0  ~ (4.1) 

2 3 

where  Kq  = 3600  cm  /cm  . Kq  was  found  to  be  independent 
of  temperature. 

(3)  The  time  dependence  of  the  volume  fraction,  at 
all  the  annealing  temperatures,  could  be  described  by  the 
empirical  equation 


Log{ 


3.08  Log{^}  - 1.679 
^0 


(4.2) 


where  t^  is  the  time  required  for  TU  recrystallization;  tg 
depends  on  the  annealing  temperature. 

(4)  The  surface  area  averaged  interface  velocity, 

Vg,  could  be  described  by  the  following  equation: 


8.5x10 


-4 


t 


■cm/  sec 


(4.3) 
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This  equation  was  obeyed  by  the  experimental  data  at  all 
the  annealing  temperatures.  Thus,  the  surface  area 
averaged  interface  velocity,  Vg,  was  found  to  be  indepen- 
dent of  temperature. 

(5)  The  activation  energy  for  the  process  was  found 
to  be  70.1  K cal/mole. 


4.2.2  Speich  and  Fisher's  Interpretation  of  Their  Data 

Speich  and  Fisher  [79]  developed  a model  to  explain 
the  experimental  behavior  of  the  surface  area  averaged 
interface  velocity,  described  by  equation  (4.3).  Their 
assumptions  are  as  follows: 

(1)  Recovery  occurs  concurrently  with  the  recrys- 
tallization process,  not  only  in  the  initial  stages,  but 
til  the  process  is  complete. 

(2)  The  decrease  in  the  stored  energy,  , due  to 
the  recovery  process  is  given  by 


dE  o 

^ = -K 
dt  r s 


(4.4) 


where  Exp{-Q^/RT}  and  is  the  activation  energy 

for  the  recovery  process,  while  is  a constant. 

(3)  The  activation  energy  for  the  recovery  process, 

, is  equal  to  the  activation  energy  for  the  growth  of 
the  recrystallized  regions. 

Using  these  assumptions,  together  with  an  expression 
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for  the  local  interface  velocity  of  the  recrystallized 
regions  developed  by  Gordon  [102],  Speich  and  Fisher  de- 
rived an  equation  similar  to  the  empirical  equation  (4.3). 
Note  that  it  is  not  possible  to  explain  the  other  experi- 
mental observations  by  using  this  model.  Furthermore, 
the  assumptions  made  by  Speich  and  Fisher,  in  their  model, 
are  physically  unrealistic. 

All  the  experimental  observations  of  Speich  and 
Fisher  will  now  be  explained  by  using  the  present  ap- 
proach. This  interpretation  does  not  require  any  of  the 
assumptions  made  by  Speich  and  Fisher. 

4 . 3 Reanalysis  of  Speich  and  Fisher’s  Data : 

Assume  that  the  recrystallized  grains  nucleate  at  the 
grain  edges  of  the  cold  worked  matrix  (Vandermeer  and 
Gordon's  work  [87]  on  the  recrystallization  in  aluminum 
supports  this  assumption)  and  they  have  a cylindrical 
shape.  Each  cylindrical  grain  reaches  a certain  length, 
£q,  soon  after  nucleation.  Once  the  grain  has  achieved 
this  length,  Zq,  its  length  does  not  increase.  Thus,  soon 
after  nucleation  a grain  grows  in  a two  dimensional 
fashion  in  the  sense  that  its  radius  increases  with  time, 
while  the  length  remains  constant.  The  radial  growth  of 
recrystallized  grains  takes  place  at  a constant  growth 
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rate,  Gj^.  It  will  be  further  assumed  that  the  nucleation 
takes  place  at  a constant  rate,  N. 

For  this  special  case,  the  parameters  in  the  model 
developed  in  section  3.2  take  the  following  values: 

P = 2;  Q = 0;  a = l 
Kv  = = Kg^  = 2tt 

{ 2 > (^-5) 

and  m=n=l;  N„  =0;  a=0 

^0 


For  these  values  of  the  parameters,  equations  (3.33), 
(3.36)  and  (3.82)  give  the  following  results: 


Vy 

(4.6) 

ex 

Sy 

2NGg  , 

(4.7) 

ex 

(Vc) 

S ex 

3£q  1 

(^2)  . 1 

(4.8) 

and 

In  order  to  compare  the  experimental  behavior  [equa- 
tions (4.1)  to  (4.3)]  with  these  predictions  of  the  model, 
it  is  necessary  to  evaluate  the  impingement  function. 

This  can  be  done  as  follows. 

For  the  two  dimensional  growth,  equation  (3.77)  gives 


V 

ex 


S 


V 

ex 


(3.77) 
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and  the  impingement  function,  is  given  by  equa- 

tion (3 . 118) : 


o dV 

«V  " -i  Sy  Exp{-^  ^ ® 


(3.118) 


Substituting  equation  (4.1)  into  the  equation  (3.118) 
gives 


(•(V  - (— 7-^)-v 


3600«n  3600i„ 

3600£q  (1 ) (1+ ) 

'V,,  [1-V^]  2 ^ g 


since  £im(i^;(V„))  0 thus  B = 0 

VI 

Furthermore,  the  function  \p(Vy)  must  be  monotonic,  and 


'P(Vy)  ->■  1 as  ->  0 . Hence, 


V 


3600£ 


1 - 


0 


= 0 


or 


£q  - 5.55  X 10  ^ cm 


(4.9) 


and 


^(V  = (1  - v^) 


(4.10) 


Substituting  equation  (4.10)  into  equations  (3.115)  and 
(3 . 116)  gives 

V, 


^ V = V 
''v 


ex 


and 


(4.11) 


- V = s 
77W  V 

(l-Vy)  ex 


(4.12) 
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Furthermore,  equations  (4.11)  and  (4.12)  give 

= (v„) 

S ^ S ^ ex 

Combining  equations  (4.6)  and  (4.11)  yields 


^ 3 


]t 


3 


(4.13) 


(4.14) 


Let  tQ  be  the  time  required  for  2%  recyrs tallization . 
Then , 


0.02 


'■  3 


]t 


3 

0 


Combining  equations  (4.14)  and  (4.15)  gives 

- 0.02(^)3 


(4.15) 


or 


Log[Tr|-]  = 3 Log[^]  - 1.699  (4.16) 

^ ’^0 

Equation  (4.16)  which  is  derived  from  the  present  model  is 
in  excellent  agreement  with  the  empirical  equation  (4.2) 
which  describes  the  experimental  behavior. 

Combining  equations  (4.8),  (4.9)  and  (4.13)  yields 

- _ 8.33  X 10'^  , 

Vg  - cm/sec  (4.17) 

Again,  equation  (4.17),  derived  from  the  model,  is  in 
excellent  agreement  with  the  empirical  equation  (4.3) 
which  describes  the  experimental  behavior. 
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Combining  equations  (2.77),  (2.11)  and  (4.12)  gives 

= 3600V^(1-V^)  (4.1) 

It  is  important  to  note  that  is  not  expected  to  change 
with  the ^annealing  temperature . The  mean  grain  intercept 
or  the  grain  size  at  the  end  of  the  recrystallization 
process  can  now  be  predicted.  Due  to  two  dimensional 
nature  of  the  growth,  the  mean  grain  intercept  of  the  re- 
crystallized grains  should  be  equal  to  (2£q) . This  result 
is  independent  of  the  number  of  recrystallized  grains  or 
their  volume  fraction.  Thus  one  can  say  that  the  grain 

size,  at  the  end  of  the  process,  should  be  equal  to  (2£q) • 
Hence , 

Recrystallized  grain  size  = 2£q  =11.1  microns (4 . 19) 

The  recrystallized  grain  size,  measured  by  Speich  and 
Fisher,  was  11  microns.  Thus  the  prediction  is  again  in 
excellent  agreement  with  the  experimental  value.  Since 
£q  is  not  expected  to  change  with  the  annealing  tempera- 
ture, the  model  predicts  that  the  recrystallized  grain 
size  should  be  independent  of  the  annealing  temperature, 
which  is  precisely  what  is  observed. 

The  present  model  explains  every  qualitative  and 
quantitative  experimental  observation.  The  assumptions 
regarding  the  constant  nucleation  rate  and  growth  rate  are 
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the  simplest  assumptions  one  can  make  and  are  supported 
by  the  theoretical  models  for  the  nucleation  and  growth  of 
the  recrystallized  regions.  The  peculiar  experimental 
behavior  of  Vg  of  the  recrystallized  grains  is  mainly 
due  to  the  two  dimensional  nature  of  the  growth.  It 
seems  that  the  complicated  model  developed  by  Speich  and 
Fisher  based  on  physically  unrealistic  assumptions  is 
rather  unnecessary. 

English  and  Backofen  [92J  have  studied  the  recrystal- 
lization behavior  in  the  same  system  and  their  data  will 
now  be  reanalyzed. 

^ ^ The  Experimental  Work  of  English  and  Backofen  [92] : 

Recyrstallization  was  studied  in  Fe-3.19%  Si  alloy. 
The  composition  of  the  alloy  was  very  close  to  the  alloy 
used  by  Speich  and  Fisher  [79].  However,  English  and 
Backofen  used  uniaxial  hot  compression  as  a mode  of  de- 
formation to  introduce  the  prestrain  rather  than  the  cold 
rolling  used  by  Speich  and  Fisher  [79].  The  volume  frac- 
tion, Vy,  and  the  surface  area,  Sy,  of  the  recrystallized 
regions  was  measured  at  different  annealing  times,  pre- 
strains and  temperatures. 

English  and  Backofen  [92]  presented  a qualitative 
and  quantitative  metallographic  evidence  for  grain  edge 
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nucleation.  Their  experimental  observations  are  as  fol- 
lows . 

(1)  The  recrystallized  grain  size  was  found  to  be 
independent  of  the  annealing  temperature. 

(2)  The  surface  area  averaged  interface  velocity, 

Vg,  could  be  described  by  the  empirical  equation 

Vg  » (4,18) 

G c 

where  c = 4.0  x 10  cm.  The  constant,  c,  was  found  to 
be  independent  of  the  annealing  temperature  but  the  con- 
stant, G,  was  temperature  dependent. 

Their  explanation  for  the  experimental  behavior  of 
Vg  is  very  similar  to  that  presented  by  Speich  and 
Fisher  [7  9]. 

4-5  Reanalysis  of  English  and  Backofen's  Data: 

It  will  be  assumed  that  the  recrystallized  grains 
nucleate  at  the  grain  edges  of  the  strained  matrix  and 
then  quickly  spread  along  the  grain  boundaries  of  the 
matrix.  The  further  growth  of  the  recrystallized  regions 
occurs  by  the  "thickening"  of  the  recrystallized  grains. 
This  amounts  to  an  assumption  of  one  dimensional  growth. 

It  will  be  further  assumed  that  all  the  nuclei  appear  at 
zero  time.  Thus  there  is  a site  saturation.  Let  Ny  be 
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the  nxamber  of  recrystallized  grains  per  unit  volume.  In 
order  to  develop  a quantitative  model,  it  will  be  assumed 
that  the  recrystallized  grains  are  of  cylindrical  shape. 
Each  grain  reaches  a certain  radius,  Rq , soon  after 
nucleation.  Once  the  grain  has  achieved  this  radius,  its 
radius  does  not  increase.  Thus,  soon  after  nucleation 
the  recrystallized  grain  grows  in  a one  dimensional 
fashion  such  that  its  length  increases  with  time,  while 
the  radius  remains  constant.  It  will  be  further  assumed 
that  the  longitudinal  growth  of  the  recrystallized  grains 
takes  place  at  a constant  rate,  G^.  Thus 

i = - t (4.19) 

where  t is  the  annealing  time  and  £ is  the  length  of  a 
recrystallized  grain. 

For  this  special  case,  the  parameters  in  the  model 
developed  in  section  3.2  take  the  following  values: 

P =0;  Q=l;  a = 0 

{N  = 0;  oj  = 3 = 1 } (4.20) 

Ky  = tt;  = Kg^  = 27t 

For  these  values  of  the  parameters,  equations  (3.33)  and 
(3.36)  give  the  following  results: 

V = (ttR^N  G )t  (4.21) 

ex  ^ ''O  ^ 


and 
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S 


V 


ex 


0 


+ (2ttRqN^  )t 


0 


(4.22) 


Combining  equations  (4.21)  and  (4.22)  gives 


S 


V 


ex 


0 0 ex 


(4.23) 


Furthermore , 


(4.24) 


To  compare  the  predictions  of  the  model  with  the  experi- 
mental behavior , it  is  necessary  to  solve  the  impingement 
problem.  This  can  be  done  as  follows. 


experimental  data.  The  behavior  can  be  satisfactorily 
described  by  a straight  line.  Thus  we  have  the  following 
empirical  equation: 


where  a and  b are  constants  and  and  are  the  volume 
fraction  and  surface  area  of  the  recrystallized  regions. 
The  values  of  a and  b for  different  values  of  prestrains 
are  given  in  Table  4.1.  Note  that  a does  not  depend  on 
the  prestrain  but  b systematically  increases  with  the  pre- 
strain. It  is  not  possible  to  verify  equation  (4.25)  for 
other  annealing  temperatures  because  English  and  Backofen 
have  not  presented  their  surface  area  data  at  other 


Figure  4.1  shows  a plot  of  (j3^)  vs.  for  the 


S 


V 


(1  - V^)(a  + b) 


(4.25) 
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V 


Fig.  4.1.  A graph  of 


vs.  Vy  at  812°C. 
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Table  4.1 

Annealing  Temperature:  812° C 

2 3 2 3 

Prestrain  b cm  /cm  a cm  /cm 


0.14 


3 


193.33 


0.30 


18 


193.33 


0.45 


21 


193.33 


annealing  temperatures.  However,  equation  (4.25)  can  be 
used  to  evaluate  the  impingement  function  at  812°C.  Com- 
bining equations  (4.25),  (4.22),  (3.115)  and  (3.116)  gives 
the  following  result  for  the  functional  formi  of  the  im- 
pingement function,  ip(V^): 


Strictly  speaking,  this  impingement  function  should  be 
used  only  for  the  annealing  temperature,  812°C  [since  only 
the  data  at  812°C  were  used  to  obtain  equation  (4.25)]. 
However , it  may  be  said  that  the  impingement  function  is 
not  expected  to  change  with  the  annealing  temperature; 
thus  equation  (4.26)  will  be  used  for  all  the  experimental 
data. 

Equations  (4.26),  (3.115)  and  (3.116)  give  the  fol- 
lowing results: 


ip(V^)  = (1  - V^)2 


(4.26) 


ex 


(4.27) 
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T " 


(4.28) 


(1-V^)  ’ex 

(Vg)  = (Vg)^^  (4.29) 

Substituting  equations  (4.27)  and  (4.28)  into  (4.23)  gives 


(1-V^)  ('^•30) 


Since  Rq  and  are  constants  for  given  temperature  and 

prestrain,  comparing  equations  (4.25)  and  (4.30)  gives 


and  b = 27tR^N„ 

u Vq 

Solving  these  equations  for  R^,  and  N„  gives 

U Vq 

R = 2 

0 (a+b) 


and 


_ b ( a+b ) ^ 


(4.32) 


(4.33) 

(4.34) 


Using  equations  (4.33)  and  (4.34)  together  with  Table  4.1 
gives  the  following  results: 


Table  4.2 

Annealing  Temperature:  812°C 


Prestrain 

Rq  (cm) 

Nt,  per  cm 
'^O 

0.14 

100x10"^ 

0.48x10^ 

0.30 

94.7x10"^ 

3.19x10^ 

0.45 

93.4x10'^ 

3.83x10^ 
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Note  that  Rq  does  not  change  appreciably  with  the  pre- 
strain but  Ny  increases  as  the  prestrain  increases. 

This  is  expected  because  the  number  of  active  nucleation 
sites  should  increase  with  increasing  strain.  The  higher 
prestrain  gives  rise  to  a finer  grain  structure  after  re- 
crystallization. However,  Rq  is  determined  by  the  grain 
structure  of  the  cold  worked  matrix  (it  is  proportional 
to  the  average  grain  face  area) ; thus  it  is  not  expected 
to  vary  with  the  prestrain  or  annealing  temperature.  In 
all  the  subsequent  calculations,  a value  of  94x10"^^  cm 
will  be  used  for  Rq . 


Combining  equations  (4.24)  and  (4.29)  and  using  R 


-4 


0 


equal  to  94x10  cm  gives 


1 (4.67x10’^) 

2g7  ^ t 


(4.35) 


This  behavior  of  Vg  is  predicted  by  the  model  is  in  ex- 
cellent agreement  with  the  experimental  behavior  described 
by  the  empirical  equation  (4.18).  For  comparison,  equa- 
tion (4.18)  is  given  as  follows: 


Thus , 


(Vg) 


1 


1 , 4x10 
G ^ T 


(4.18) 


G = 


(4.36) 
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Since  depends  on  the  annealing  temperature,  G should 
also  depends  on  the  annealing  temperature.  This  is  pre- 
cisely what  is  observed.  On  the  other  hand,  the 


coefficient  of  (1/t)  in  the  equation  (4.35)  is  (-^)  and 
hence  it  is  not  expected  to  change  with  the  prestrain  or 
annealing  temperature.  Again,  this  is  in  excellent  agree- 
ment with  the  observed  behavior. 

Combining  equations  (4.21)  and  (4.27)  gives 

V, 


V 


V 


[ttRoNv^Gl]  t 


(4.37) 


Substituting  equation  (4.32)  into  (4.37)  gives 

V, 


V 


V 


[2bGj^]  t 


(4.38) 


Since  the  quantity  in  the  square  bracket  is  independent 
of  time,  the  model  predicts  that  the  plot  of  (-■ — ^)  vs.  t 
should  be  linear  and  the  straight  line  should  pass  through 
the  origin.  Figure  4.2  shows  a plot  of  ^ vs.  t for 
the  experimental  data  at  812°C.  The  plot  is  linear  and 
the  straight  line  goes  through  the  origin.  Thus  the  ex- 
perimental behavior  is  in  excellent  agreement  with  the 
predictions  of  the  model.  Let  t^  be  the  time  required 
for  total  recrystallization.  Then,  equation  (4.38)  gives 

0. 11  = [2bGj^]tQ  (4.39) 

Substituting  equation  (4.39)  into  (4.38)  gives 


0.45 
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Figure 
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TT^=0-11[^]  (4.40) 

V ^0 

V 

Thus,  the  model  predicts  that  a plot  of  vs.  (— ) 

should  be  linear  and  the  data  points  at  all  the  prestrains 
and  annealing  temperatures  should  fall  on  a single 
straight  line  of  slope  0.111.  Furthermore,  the  straight 
line  should  pass  through  the  origin.  Figure  4.3  shows 
such  a plot  for  the  experimental  data  at  various  pre- 
strains and  annealing  temperatures.  Again  the  agreement 
between  the  experimental  behavior  and  prediction  of  the 
model  is  very  good.  Equation  (4.39)  can  now  be  used  to 
evaluate  the  values  of  the  growth  rate  constant,  , if  b 
is  known.  Table  4.1  gives  values  of  b at  three  pre- 
strains. It  may  be  assumed  that  b is  independent  of 
temperature.  The  values  of  evaluated  in  this  manner 
are  given  in  the  Table  4.3.  The  first  three  values 

of  G^  in  Table  4.3  can  be  used  to  evaluate  the  activation 
energy  for  the  growth  process  (this  is  because  the  pre- 
strain is  not  very  different  for  these  data) . A plot  of 
In  Gj^  vs.  1/T  is  shown  in  Figure  4.4.  The  plot  is  linear, 
indicating  an  Arrhenius  nature  of  the  temperature  de- 
pendence of  Gj^.  The  activation  energy  may  be  evaluated 
from  the  slope  of  this  plot  and  the  calculated  value  is 
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Figure  4.3  A graph  of  Vy/Cl-Vv)  vs.  t/tQ  at  different 
temperatures  and  prestrains. 
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Per  ® k 


Figure  4.4.  An  Arrhenius  plot  of  £n  vs . 1/T. 
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Table  4.3 


Prestrain 

Annealing  Temperature  °C 

cm/ sec 

0.39* 

750 

0.230x10"^ 

0.45 

812 

1.32x10'^ 

0.39* 

911 

5.87x10"^ 

0.26** 

750 

0.086x10“^ 

0.30 

812 

0.56x10"'^ 

*b  value  assumed  to  be  same  as  that  for  the  pre- 
strain of  0.45 

**b  value  assumed  to  be  same  as  that  for  the  pre- 
strain of  0.30 

as  follows : 

rActivation  energy.  = Qp  = 47  K cal/mole 
for  growth  ^ 

This  simple  model  explains  all  the  experimental  be- 
havior quantitatively.  The  analysis  clearly  demonstrates 
the  applicability  of  the  present  approach.  A similar 
experimental  behavior  has  been  observed  in  some  other  re- 
crystallization studies  [93,94]  and  can  be  explained  by 
using  the  present  approach.  The  next  two  chapters  deal 
with  the  experimental  studies  done  during  the  present 
work,  and  it  will  again  be  shown  that  this  approach  is 
very  successful  in  analyzing  the  experimental  data. 


CHAPTER  V 

NUCLEATION  AND  GROWTH  OF  AUSTENITE  FROM  PEARLITE 
5 . 1 Introduction : 

The  formation  of  austenite  has  been  studied  in  much 
less  detail  than  its  decomposition.  The  previous  investi- 
gations are  not  entirely  satisfying.  In  some  of  the 
studies,  quantitative  data  is  not  available,  while  in 
other  cases  there  are  some  contradictions. 

This  chapter  deals  with  an  experimental  investigation 
of  the  formation  of  austenite  from  pearlite  in  a eutectoid 
Fe-C  alloy.  A brief  survey  of  the  previous  work  is  pre- 
sented in  the  first  section  and  the  subsequent  sections 
deal  with  the  experimental  work  and  data  analysis. 

5 . 2 Literature  Survey: 

A study  of  the  austenitization  of  pearlite  was 
carried  out  by  Roberts  and  Mehl  [105]  in  eutectoid  steels. 
They  concluded  that  the  overall  transformation  could  be 
described  satisfactorily  by  a constant  nucleation  and 
growth  rate  model.  Speich  and  Szirmae  [106]  studied 
formation  of  austenite  from  pearlite  in  hypereutectoid 
Fe-C  alloys.  They  claim  that  the  nucleation  of  austenite 
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occurs  preferentially  at  the  intersections  of  the  pearlite 
colonies.  The  following  observations  were  made  by  Speich 
and  Szirmae  [106]. 

(1)  The  relationship  between  the  volume  fraction  of 
austenite,  Vy,  and  its  surface  area,  Sy,  could  be  de- 
scribed by  an  empirical  equation: 

sj  - vj[l-v}l  - A2{vJ(1-vJ)}2  (5.1) 

where  and  A2  depend  on  temperature  but  not  on  the 
volume  fraction. 

(2)  The  surface  area  averaged  interface  velocity  was 
found  to  be  constant  at  a given  austenitizing  temperature. 

(3)  The  time  dependence  of  the  volume  fraction 
could  be  described  by  the  following  equation: 

= 1 - Exp{-B^t^}  (5.2) 

where  is  a kinetic  constant.  From  these  three  observa- 
tions Speich  and  Szirmae  concluded  that  the  transformation 
takes  place  by  the  site  saturation  of  the  austenite  nuclei 
along  intersections  of  the  pearlite  colonies  and  the  sub- 
sequent three  dimensional  growth  of  the  nuclei  at  a con- 
stant rate. 

Note  that  equation  (5.2)  assumes  a random  impinge- 
ment. For  such  a case  of  site  saturated  three  dimensional 
growth,  the  model  developed  in  section  3.2  gives 
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(1-Vj) 


ex 


(5.3) 


where  is  a kinetic  constant.  Thus  a plot  of 

sr 

vs.  £n  t should  be  linear  with  a slope  of  2 0 

(i-vj) 

Such  a plot  for  Speich  and  Szirmae's  data  [106]  was  found 
to  be  linear  but  the  slope  was  observed  to  be  0.90  instead 
of  2.0.  Furthermore,  the  surface  area  averaged  interface 
velocity,  Vg , calculated  from  equations  (5.1)  and  (5.2) 
was  found  to  increase  systematically  with  time,  increasing 
almost  by  one  order  of  magnitude.  Thus  the  conclusions 
reached  by  Speich  and  Szirmae  [106]  need  to  be  re- 
examined . 

Karlson  [107]  has  studied  the  transformation  kinetics 
of  the  austenitization  process  in  Fe-0.18%C  alloy  at 
855  C.  It  was  observed  that  pearlite  transforms  to 
austenite  in  less  than  one  second  and  the  further  growth 
of  austenite  occurs  by  consumption  of  ferrite  by  the 
existing  austenite  grains.  Karlson  [107]  concluded  that 
diffusion  of  carbon  through  austenite  is  the  rate  con- 
trolling step. 

The  formation  of  austenite  has  been  studied  by 
several  other  workers  [108-112]  but  most  of  these  papers 
deal  with  either  the  effect  of  rapid  heating  on  the  forma- 
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tion  of  austenite  or  the  experiments  are  done  on  steels 
with  significant  amounts  of  alloying  elements.  A review 
article  on  the  austenitization  process  is  given  by 
Paxton  [113] . 


5 . 3 Experimental  Work: 

The  experiments  were  done  on  Fe-0.837cC  alloy  supplied 
by  Bethlehem  Steel  Corporation.  The  alloy  composition  is 
given  in  Table  5.1. 


Table  5.1 


Element 

Alloy  Composition 
Wt.  Percent  Element 

Wt . Percent 

C 

0.83 

Mo 

0.04 

Mn 

<0.01 

Cu 

0.001 

Si 

<0.002 

Sn 

<0.002 

S 

<0.01 

A1 

<0.015 

Cr 

0.01 

Ti 

<0.001 

Ni 

<0.01 

N 

0.001 

The  alloy  was  homogenized  at  1000°C  for  three  days. 
Zirconium  chips  were  used  to  minimize  the  oxidation  and 
decarburization  during  the  homogenizing  treatment.  Three 
samples  were  cut  from  the  different  portions  of  the  homo- 
genized bars  and  were  polished  for  the  microscopic  exam- 
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All  tli6  samples  were  found  to  be  completely 
pearitic  and  there  were  no  concentration  gradients. 

The  homogenized  bar  was  machined  from  to  %"  dia- 
meter and  the  samples  of  1/16”  thickness  were  cut  from 
the  diameter  rod.  These  samples  were  used  for  the 
experimental  work. 

5.3.1  Starting  Structure  for  the  Austenitization  Studies 
A uniform  inter lamellar  spacing  was  produced  in  the 
samples  by  the  following  heat  treatment. 

The  specimens  were  austenitized  at  850°C  for  fifteen 
ininutes  in  a salt  bath  furnace  and  then  quickly  trans- 
ferred to  another  salt  bath  at  680°C.  The  samples  were 
kept  at  this  temperature  for  twenty  minutes  and  then 
quenched  in  iced  brine.  The  salt  baths  were  sprinkled 
with  graphite  powder  to  avoid  decarburization. 

Three  samples  were  polished  and  nital  etched.  The 
structure  of  these  samples  was  found  to  be  completely 
pearlitic.  The  pearlite  spacing  was  observed  to  be  uni- 
form and  the  same  in  all  three  samples. 

5.3.2  Austenitization  Experiments 

The  specimens  were  austenitized  at  750°C  for  various 
times  in  a salt  bath  furnace  and  then  quenched  in  iced 
brine.  The  temperature  control  of  the  furnace  was 
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approximately  ±2°C.  The  procedure  was  repeated  at  730°C 
and  740°C,  The  experiments  had  to  be  repeated  several 
times  to  get  well-spaced  austenitization  times  and  uniform 
microstructures . 


5 . 4 Metallography : 

Each  sample  was  cut  in  the  middle  and  mounted  "edge 
on."  This  "edge  on"  mounting  was  done  to  avoid  any  sur- 
face effects.  The  polishing  was  done  by  the  standard 
polishing  and  grinding  techniques.  The  final  wheel 
polishing  was  done  by  using  6 micron  diamond  paste. 

Before  etching,  each  sample  was  examined  under  the  micro- 
scope. The  relief  between  pearlite  and  martensite  was 
small.  The  samples  were  etched  with  a 2%  nital  solution. 
The  representative  microstructures  are  shown  in  Fig- 
ures 5.1  to  5.3. 

It  was  observed  that  the  austenite  grains  form  uni- 
formly in  the  pearlite  matrix  and  at  low  volume  fractions 
the  microstructure  on  the  section  is  typically  one  of  the 
polydispersed  spheres  distributed  in  the  matrix.  ,Some 
J^Fl§H£ite  was  left  even  after  the  completion  of_the_ 
austenitization  process.  In  most  of  the  cases,  this 
Xesidyal^cementite  was  in  the  form  of  spherodite  rather 
than _lame liar  shape.  Thus  one  may  conclude  that  the  com- 
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Figure  5.1.  A photomicrograph  of  a sample  partially 
austenitized  at  750°C  for  23  sec. 
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Figure  5.2.  A photomicrograph  of  a sample  partially 
austenitized  at  750°C  for  32  sec. 
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Figure  5.3.  A photomicrograph  of  a sample  partially 
austenitized  at  750°C  for  45  sec. 
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position  of  austenite  at  the  end  of  the  austenitization 
process  is  not  uniform. 

The  volume  fraction  and  surface  area  of  martentsite 
(which  was  austenite  at  the  reaction  temperature)  was 
measured  in  each  sample  by  using  the  standard  stereolo- 
gical  techniques.  The  measurements  were  carried  out 
manually.  The  image  was  projected  on  the  T.V.  screen  of 
the  Quantimet  720  and  the  measurements  were  performed  by 
superimposing  a grid.  At  least  fifty  grid  placements  were 
used  in  each  sample.  The  volume  fraction  was  measured 
by  point  counting  while  the  surface  area  was  measured 
by  the  line  intercept  count.  The  measured  values  of  the 
f^3-Ction  and  surface  area  of  austenite  are  given 
in  Tables  5.2  to  5.4. 


5 . 5 Data  Processing: 

In  this  case,  one  may  assume  that  the  impingement  is 
random.  Thus,  it  is  possible  to  evaluate  the  extended 
properties  by  using  equations  (2.15)  and  (2.16).  The  ex- 
tended properties  evaluated  in  this  manner  are  given  in 
Tables  5.2  to  5.4. 

A plot  of  Sy  vs.  Vy  is  shown  in  Figure  5.4.  The 
ex  ex 

plot  is  nonlinear.  Such  a plot  is  expected  to  be  linear 
for  a two  or  one  dimensional  growth  [see  equations  (3.77) 


Table  5.2 

Global  Properties  of  Austenite 
stenitizing  Temperature:  750° 
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Table  5.3 

Global  Properties  of  Austenite 
■Stenitizing  Temperature;  740° 


93 


CJ 


p 

< 


B 

o 


B 

o 


>-  > 


> 

> 


ro 


6 

o 


CSJ 

B 

a 

> 


>-  > 
> 


60 

C 

•P^ 

• 

N 

O 

*P^ 

Q) 

4J 

W 

•H 

'w' 

c 

01 

0) 

B 

w 

•H 

p 

•U 

<C 

o 


CM 

o 

CO 

00 

o 

o 

VO 

ON 

CM 

CO 

CM 

00 

CO 

o^ 

Ov 

1— 1 

CM 

O 

O 

O 

T— 1 

r-- 

CO 

VO 

•<}■ 

I— 1 

CM 

o 

o 

o 

o 

o 

t— 1 

CM 

VO 

CM 

1—1 

oo 

VO 

CO 

rH 

00 

00 

CM 

UO 

VO 

LiO 

CO 

O 

vX> 

OO 

CO 

o 

rH 

CO 

UO 

00 

o 

O 

O 

o 

o 

o 

o 


o 

m 


o 

VO 


o 


o 

00 


o 

Ov 


CN  CO 


m VO 


Table  5.4 

Global  Properties  of  Austenite 
Austenitizing  Temperature:  730° 
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Figure  5.4.  A graph  of  vs.  V^0jj  during  austenitiza- 

tion of  pearlite. 
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t 

Figure  5.5.  A graph  of  fS^dt  vs. 


at  750°C. 
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S/10  -t 

J * 

0 


sec./cm. 


t 

Figure  5.6.  A graph  of  /S„dt  vs.  V„ 

0 V V 


at  740°C. 
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t 

Figure  5.7.  A graph  of  /sXdt  vs.  at 
730°C.  o''  '' 
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and  (3.90)].  Thus,  it  may  be  concluded  that  the  austenite 

growth  occurs  in  a three  dimensional  manner. 

t 

A plot  of  {/Sydt}  vs.  Vy  was  found  to  be  linear  at 
all  three  austenitizing  temperatures.  These  plots  are 
shown  in  Figures  5.5  to  5.7.  The  slope  of  the  straight 
line  changes  with  the  austenitizing  temperature.  Thus, 
we  have  the  following  empirical  equation: 

t 

/Sjdt  - a(T)  • (5.4) 

where  a(T)  does  not  depend  on  the  time  but  it  decreases 

as  the  austenitizing  temperature  increases . 

2 

A plot  of  {S^  /t  } vs.  t is  shown  in  Figures  5.8  to 
ex 

5.10.  This  plot  is  also  linear  at  all  the  austenitizing 
temperatures.  Thus,  the  following  empirical  equation  de- 
scribes the  time  dependence  of  the  extended  surface 
area. 

Sj  = t^[B|  + B^t]  (5.5) 

ex 

where  B^  and  B^  depend  on  the  temperature.  Table  5.5  gives 
the  values  of  a(T),  B^  and  B^  at  different  austenitizing 
temperatures . 

5 . 6 Data  Analysis : 

The  nonlinear  plot  of  vs.  (see  Figure  5.4) 

'^ex  ^ex 
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Y 2 

A graph  of  SA  /t  vs . t at 
750°C.  ^ex 


Figure  5.8. 


101 


Figure  5.9.  A graph  of  sj  /t^  vs.  t at  740°C. 

''ex 


<-< 
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Figure  5.10.  A graph  of  sl 
730°C.  ^ 


2 

/t  vs.  t at 
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Table  5.5 


Austenitizing 
Temperature  °C 

a(T)  sec/cm 

Y 2 

1/cm-sec 

B2  1/cm-sec^ 

730 

1.248x10^ 

0.102x10"^ 

0.036x10'^ 

740 

3.2x10^ 

1.35x10'^ 

4.023x10"^ 

750 

1.83x10^ 

0.0528 

0.026 

indicates  that  the  austenite  growth  is  three  dimensional. 
Since  the  matrix  (pearlite)  composition  is  not  expected 
to  change  during  the  process,  it  can  be  said  that  the 
austenite  growth  rate  is  only  size  dependent-.  Thus,  equa- 
tion (3.47)  can  be  used  to  deduce  the  mathematical  form  of 
the  austenite  growth  rate.  This  can  be  done  as  follows. 
Equation  (5.4)  can  be  rewritten  in  the  following 

form: 

t t dVj 

/sjdt  = a(T)  / ^ dt  (5.6) 

0 0 

Substituting  equations  (2.15)  and  (2.16)  into  this  equa- 
tion yields 

t dVj 

/ e - a(T)  _^}dt  = 0 (5.7) 

Equation  (5.7)  can  be  satisfied  for  any  arbitrary  value 
of  t only  if  the  integrand  is  identically  equal  to  zero. 
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Thus  , 


dV 


a(T) 


V 


ex 


ex 


dt 


= 0 


or 


fS^  dt  = a(T)  Vj 
0 ex  ex 


(5.8) 


Now,  equation  (3.47)  gives 


^ V t 

/V(t-T)  (x)dT  = /S(t-T)  Vj  (x)dT 


0 ^ex  0 ^ex 


(3.74) 


Integrating  the  left-hand  side  of  this  equation  by  parts 
yields 


^^_9V(^t^-x)^x  . /sj  (t)dx  dx  = /S(t-x)V:j:  (x)dx(5.9) 

U 0 ex  0 


ex 


Substituting  equation  (5.8)  into  (5.9)  gives 
t 


■Y 


0 ex 


(x)[S(t-x)  - a(T)  (-^^^^-l^)^]dx  = 0 (5.10) 


Since  the  expression  in  the  square  bracket  is  a function 
of  (t-x)  only,  the  integrand  has  to  be  identically  equal 
to  zero  if  equation  (5.10)  is  to  be  satisfied  at  any  value 
of  time,  t.  Thus 

S(t-x)  - a(T)  (H^lll)^  = 0 
or 

1 .9V(t-x).,  ^ 1 

S(t-x)  ••  at  Jx 


(5.11) 
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The  expression  on  the  left-hand  side  of  equation  (5.11) 
is  the  surface  area  averaged  interface  velocity  of  one 
P^^ticle  (nucleated  at  time  x).  The  expression  is  valid 
for  any  value  of  x ; thus  it  is  valid  for  all  the  parti- 
cles . 

It  has  already  been  deduced  that  the  growth  is  three 
dimensional.  It  can  be  qualitatively  deduced  from  the 
microstructures  that  the  austenite  particles  have  a 
spherical  shape.  Let  R be  the  radius  of  the  austenite 
particle  at  time  t whose  nucleation  time  is  x.  One  can 
write 


V(t-x)  = ^ 


and 


(5.12) 


Equations 


Combining 


S(t-x)  = 4tt 

(5.12)  and  (5.13)  give 

, 1 , r3V(t-x).  ^ 

S(t-x)  8t  -*x 

equations  (5.11)  and  (5.14 


[—1  = 

*-dt-*x  a (T) 


dR, 

dt^x 

yields 


(5.13) 


(5.14) 


(5.15) 


Thus  the  austenite  growth  occurs  at  a constant  growth  rate 
and  its  value  at  a particular  temperature  is  given  by 
(l/a(T)).  Thus,  we  have  deduced  the  growth  rate  of 


austenite  from  only  the  time  variation  of  the  volume  frac- 
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tion  and  surface  area  without  any  knowledge  or  assumption 
regarding  the  nucleation  rate.  To  the  best  knowledge  of 
the  author,  this  has  never  been  done  before.  The  result 
clearly  demonstrates  the  potential  and  applicability  of 
equations  (3.47)  to  (3.49)  developed  in  the  third 
chapter . 

The  nucleation  rate,  N(t) , will  now  be  evaluated. 

For  a particle  nucleated  at  t = 0 , equations  (5.13)  and 
(5.15)  give  the  following  result; 

S(t)  = • t^  (5.16) 

[a(T)]^ 

Taking  the  Laplace  transform  of  both  the  sides  of  this 
equation  yields 

L{S(t)}  = (5.17) 

[amr 

Note  that  S on  the  right-hand  side  of  this  equation  is  the 
parameter  of  the  Laplace  transform  and  it  should  not  be 
confused  with  the  surface  area. 

Similarly,  equation  (5.5)  gives 


2b{  6B^ 

L{S^  (t)}=-^+-^  (5.18) 

ex  S S 

Combining  equations  (3.48),  (5.17)  and  (5.18)  yields 


N(t)  = L'^{2b][  + 


6B): 


(5.19) 
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or 

2 

N(t)  = J {B]^6(t)  + 3Bp  (5.20) 

where  (t)  is  the  delta  function  and  it  is  identically 
equal  to  zero  for  t 0.  Thus  for  t 0,  one  can  write 

N(t)  = ^ bJ  (5.21) 

Thus,  the  nucleation  rate  is  constant  and  its  value  is 
given  by  equation  (5.21).  The  values  of  the  nucleation 
rate  and  growth  rate,  evaluated  by  using  equations  (5.21) 
and  (5.15),  are  given  in  Table  5.6 

Table  5.6 

Temperature  °C  Growth  rate  cm/ sec  Nucleation  rate  1/cm^-sec 


730 

0.80x10"^ 

1.34x10 

740 

0.31x10'^ 

9.8x10^ 

750 

0.55x10"^ 

20x10^ 

Brandt  [114]  has  solved  the  diffusion  equations  for 
the  growth  of  pearlite,  assuming  that  the  volume  diffusion 
of  carbon  is  the  rate  controlling  process.  For  the  alloy 
of  eutectoid  composition,  similar  technique  can  be  used 
to  evaluate  the  growth  rate  of  austenite  from  pearlite. 
This  method  was  used  by  Speich  and  Richards  [116]  to 
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evaluate  the  growth  rate  of  austenite  between  820°C  and 
900°C.  Evaluation  of  the  growth  rate  by  using  this  tech- 
nique requires  a solution  to  four  nonlinear  simultaneous 
equations  and  it  gives  a complicated  expression  for  the 
temperature  dependence  of  the  austenite  growth  rate. 
However,  at  a given  temperature  a constant  growth  rate  is 
predicted,  which  is  in  agreement  with  the  present  experi- 
mental observations. 


5 . 7 Conclusions : 

(1)  The  growth  rate  of  austenite  was  deduced  from 
the  time  variation  of  its  volume  fraction  and  surface 
area,  without  any  knowledge  or  assumptions  regarding  the 
nucleation  rate.  The  result  clearly  demonstrates  the 
potential  and  the  practical  feasibility  of  the  present 
approach. 

(2)  The  austenite  growth  rate  was  found  to  be  con- 
stant at  a given  temperature. 

(3)  The  nucleation  rate  of  austenite  was  evaluated 
from  the  time  dependence  of  its  extended  surface  area. 
The  nucleation  rate  was  found  to  be  constant  at  a given 
temperature . 


CHAPTER  VI 

FORIIATION  OF  AUSTENITE  FROM  "SPHERODIZED”  CEMENTITE 
AND  FERRITE  AGGREGATES 

6 . 1 Introduction : 

There  are  very  few  experimental  studies  on  the  forma- 
tion of  austenite  from  spherodized  cementite  structures. 
These  previous  investigations  are  incomplete  and  incon- 
clusive. This  chapter  deals  with  the  experimental  study 
of  the  formation  of  austenite  from  spherodized  cementite 
structures,  in  a eutectoid  Fe-C  alloy.  The  investigation 
clearly  demonstrates  the  effect  of  "geometry"  on  the 
kinetic  and  kinematics  of  the  phase  transformation.  The 
experiments  were  carried  out  on  three  different  initial 
distributions  of  spherodized  cementite  and  at  three  dif- 
ferent austenitizing  temperatures.  The  data  are  analyzed 
in  a systematic  manner  (which  is  the  most  serious  weak- 
ness of  the  previous  studies)  by  using  the  approach  de- 
veloped in  the  previous  chapters.  A brief  literature 
survey  is  presented  in  the  first  section.  The  second 
section  deals  with  the  experimental  work;  the  subsequent 
sections  are  concerned  with  the  data  processing  and 
analysis . 
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6 . 2 Literature  Survey: 

Roberts  and  Mehl  [105]  studied  nucleation  and  growth 
of  austenite  from  spherodized  cementite  structures.  They 
observed  that  the  nucleation  of  austenite  does  not  occur 
at  every  carbide  particle  because  a single  austenite 
grain  may  envelope  several  carbide  particles.  Roberts 
and  Mehl  [105]  observed  that  the  ferrite  grain  boundaries 
play  a very  important  role  in  the  austenitization  of  the 
spherodized  cementite  structures.  A qualitative  de- 
scription of  the  austenitization  process  is  given  by 
Bain  [117].  Molinder  [118]  studied  austenitization  of  a 
spherodized  cementite  structure  in  1.27  wt.%  carbon 
steel.  Molinder  [118]  concluded  that  the  process  is 
interface  controlled.  The  kinetics  of  austenitization 
process  was  studied  by  Neizweidz,  Partom,  Taub  and 
Weiss  [119]  in  a graphitic  steel.  The  transformation 
was  studied  by  measuring  the  volume  fraction  of  austenite 
and  cementite  at  various  times  during  the  process.  The 
data,  on  volume  fraction  of  austenite,  were  analyzed  by 
using  an  equation  of  the  type 

= 1 - Exp{-Bt’^}  (6.1) 

The  values  of  n varied  from  0.8  to  1.74  during  different 
stages  of  the  transformation  at  a given  temperature. 
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Their  volume  fraction  data  were  reanalyzed  during  the 
present  work  and  it  was  found  that  following  simple  equa- 
tion describes  the  time  variation  of  the  volume  fraction 
of  austenite: 


-Bt 

(i-v^)  ®i^ 


(6.2) 


where  is  a function  of  temperature  only.  It  was  shown 
in  the  previous  chapters  that  if  the  impingement  function 


is  described  by  (1-V^)'^,  then 


V. 


V. 


V 


^ex 


Combining  equations  (6.2)  and  (6.3)  gives 


V 


V ®1^ 
ex 


(6.3) 


(6.4) 


Thus,  under  the  assumption  of  equation  (6.3),  equa- 
tion (6.4)  describes  a linear  relationship  between  the 
extended  volume  fraction  of  austenite  and  time.  Equa- 
tion (6.4)  may  be  interpreted  as  the  site  saturated  two 
dimensional  parabolic  growth  of  austenite.  Thus,  the 
complicated  numerical  solution  to  the  diffusion  problem, 
presented  by  Neidzweidz  et  al . [119] , may  not  be 

necessary . 

Judd  and  Paxton  [120]  have  studied  austenitization 
of  spherodized  cementite  and  ferrite  aggregates  in  some 
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low  carbon  steels.  They  have  developed  a mathematical 
model  for  the  kinetics  of  dissolution  of  spherodized 
cementite.  The  model  does  not  take  into  account  the  geo- 
metrical impingement  between  the  growing  austenite 
grains;  thus,  it  is  applicable  only  if  the  volxome  frac- 
tion of  spherodized  cementite  is  small  and  the  inter- 
particle distance  is  large.  Furthermore,  the  model  as- 
sumes that  the  spherodized  cementite  is  in  the  form  of 
discrete  particles.  It  will  be  demonstrated  that  this 
may  not  be  necessarily  true  if  the  cementite  is  situated 
at  the  ferrite  grain  edges.  Judd  and  Paxton  [120]  con- 
cluded that  the  volume  diffusion  of  carbon  in  austenite 
is  the  rate  controlling  step  for  the  process . 

Speich  and  Szirmae  [106]  studied  the  formation  of 
austenite  from  spherodized  cementite  structure  in  some 
Fe-C  alloys . They  observed  that  growth  of  austenite 
grains  is  inhibited  across  the  ferrite  grain  boundaries. 
This  may  suggest  an  orientation-relationship  between  the 
austenite  and  ferrite  grains.  A review  article  on  the 
austenitization  process  is  given  by  Paxton  [113]. 

Mathematical  models  for  dissolution  of  spherodized 
cementite,  during  the  austenitization  process,  have  been 
developed  by  Neidzweidz,  Weiss,  Partom  and  Taub  [119]  and 
Judd  and  Paxton  [120].  Molinder  [118]  has  developed  a model 
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for  the  time  variation  of  the  average  austenite  composi- 
tion during  the  process. 

6 . 3 Experimental  Work: 

The  experiments  were  done  on  the  Fe-0.83  wt.%  C 
alloy  supplied  by  Bethlehem  Steel  Corporation.  The  com- 
position of  the  alloy  is  given  in  Table  5.1.  The  homo- 
genization treatment  is  described  in  section  5.2.  The 
samples  of  diameter  and  1/16”  thickness  were  used  for 
the  experiments. 

6.3.1  Starting  Structure 

The  samples  were  austenitized  at  850°C  for  thirty 
minutes  in  a salt  bath  furnace  and  then  quenched  in  the 
iced  brine.  The  microstructure  of  three  samples  was  ex- 
amined and  was  found  to  be  completely  martensitic.  The 
hardness  of  each  quenched  sample  was  measured  on  the 
Rockwell  ”c”  scale  and  the  samples  with  the  average  hard- 
ness less  than  Rc.  63  were  discarded  (the  hardness  of 
eutectoid  martensite  is  Rc . 65).  The  martensitic  samples 
were  vacuum  sealed  in  vycor  tubes  and  annealed  in  a muffle 
furnace  at  690°C  to  produce  the  spherodized  cementite 
structures.  Three  different  initial  distributions  of 
spherodized  cementite  were  produced  by  annealing  the 
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samples  for  different  annealing  times.  The  details  of 
annealing  times  and  the  letter  designations  for  these 
starting  structures  are  given  in  Table  6.1. 

Table  6.1 

Annealing  Temperature:  690°C 

Annealing  Time  Letter  Designation 
1 day 
1 week 
1 month 

6.3.2  Austenitizing  Experiments 

The  samples  of  all  the  three  batches,  D^,  D2  and 
D^,  were  partially  austenitized  at  750°C  for  various 
times  in  a salt  bath  furnace  and  then  quenched  in  iced 
brine.  The  temperature  control  of  the  surface  was  ±2°C. 
The  remaining  samples  of  batch  D2  were  partially 
austenitized  at  740°C  and  760°C  for  various  times  and  then 
quenched  in  iced  brine. 

6.3.3  Metallography 

The  polishing  was  done  by  standard  grinding  and 
polishing  techniques.  The  samples  were  ground  through  a 
series  of  SiC  papers  to  600  grit  and  then  polished  on  a 
polishing  wheel  with  six  micron  diamond  paste.  In  order 
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to  measure  the  global  properties  of  cementite  phase  by- 
using  the  Quantimet  720,  the  samples  -were  etched  by  using 
the  etching  procedure  [121]  described  in  Table  6.2. 

Table  6.2 

Etching  Procedure  for  Etching  Spherodized  Cementite 
Composition  of  etching  solution: 
picric  acid:  2 gm 

Sodium  hydroxide:  24.5  gm 

Water:  73.5  cc 

Etch  the  sample  for  ten  to  fifteen  minutes  in  the 
boiling  etching  solution 

In  the  etched  structures  spherodized  cementite  was  colored 
dark  brown,  while  ferrite  and  martensite  were  unaffected. 
This  etching  solution  does  not  etch  the  ferrite  grain 
boundaries.  A typical  microstructure  revealed  in  this 
manner  is  shown  in  Figure  6.1.  The  volume  fraction, 
vj,  surface  area,  S^,  and  total  curvature,  of  the 

spherodized  cementite  was  measured  in  all  the  samples  by 
using  the  image  analyzing  computer,  Quantimet  720.  De- 
tails of  the  procedure  is  given  elsewhere  [56] . The 
values  of  the  global  properties  of  cementite  are  given  in 
Tables  6.3  to  6.8.  In  addition  to  these  global  prop- 
erties, the  inverse  cumulative  intercept  distribution  of 
the  linear  intercepts  of  cementite  phase  was  measured  in 
each  sample.  This  is  explained  as  follows. 


116 


Figure  6.1.  A photomicrograph  of  the  starting  structure 
D3  etched  with  the  boiling  sodium  picrate 
solution . 
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Global  Properties  of  Cementite  During  Austenitization  at  750 
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Table  6.7 

Global  Properties  of  Cementite  During  Austenitization  at  760 

Starting  Structure: 
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Table  6.8 

Global  Properties  of  Cementite  During  Austenitization  at  740 

Starting  Structure:  D., 
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Let  nj^(£)  be  the  distribution  function  of  the  linear 
intercepts  such  that  nj^(£)d£  gives  the  number  of  inter- 
cepts per  unit  length  of  the  test  line,  whose  length  is 
between  £ and  (£  + d£) . The  inverse  cumulative  distri- 
bution function  of  the  linear  intercepts,  Nj^^(£),is  de- 
fined by  the  following  equation: 

m 

N^^(£)  = /nj^(£)d£  (6.5) 

where  £^  is  the  maximum  intercept  length  observed.  Thus, 
Nl>(^-)  is  the  number  of  intercepts  of  length  greater  than 
£,  per  unit  length  of  the  test  line.  It  is  possible  to 
measure  N^^(£)  for  different  values  of  £ by  using  the 
Quantimet  720.  The  Nj^^(£)  function  for  spherodized 
cementite  was  measured  in  each  sample.  The  values  of 

for  different  values  of  £ are  given  in  Tables  6.9 
to  6.14.  If  the  particles  are  discrete  and  have  a simple 
convex  shape,  then  it  is  possible  to  evaluate  the  three 
dimensional  size  distribution  from  the  measured 
Nl>(0  [51]. 

One  starting  structure  (batch  D2)  was  polished  and 
etched  with  27o  nital  and  the  volume  fraction  and  surface 
area  of  the  spherodized  cementite  was  measured  manually. 
The  sample  was  repolished  and  the  structure  was  revealed 
by  the  etching  procedure  described  in  Table  6.2.  The 
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volume  fraction  and  surface  area  of  the  spherodized 
cementite  were  measured  by  using  the  Quantimet.  The 
values  of  the  volume  fraction  and  surface  area  of  the 
spherodized  cementite  measured  by  these  two  procedures 
are  given  in  Table  6.15. 

Table  6.15 

Starting  Structure:  D2 

2%  nital  etch  and  Boiling  sodium  picrate  etch  and 
manual  measurements Quantimet  measurements 

Sy  cm  / cm  Vy  Sy  cm  /cm 

0.128  1456  0.13  1500 

No  significant  difference  is  observed  in  the  values  of 
the  volume  fraction  and  surface  area  measured  on  the 
same  sample,  etched  by  two  different  etching  procedures. 
Thus , one  may  conclude  that  the  etching  technique 
(Table  6.2)  did  not  give  rise  to  any  artifacts. 

It  is  not  possible  to  reveal  the  ferrite  grain 
boundaries  by  using  the  etching  procedure  described  in 
Table  6.2  In  order  to  etch  both  the  spherodized 
cementite  and  the  ferrite  grain  boundaries,  the  fol- 
lowing etching  procedure  was  used.  The  samples  were 
first  etched  by  using  the  procedure  described  in  Table  6.2 
and  then  re-etched  by  using  the  etching  procedure  de- 
scribed in  Table  6.16. 
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Table  6.16 

Etching  Procedure  for  Etching  the  Ferrite 
Grain  Boundaries  and  Martensite 

Solution  A:  2 gm  picric  acid  in  100  cc  of 

methyl  alcohol 

Solution  B:  Tk  nital 

Immerse  the  sample  in  Solution  A for  5 seconds  and  then 
immerse  it  in  Solution  B for  5 to  10  seconds. 

A typical  microstructure  of  the  starting  as-coarsened 
samples  revealed  in  this  manner  is  shown  in  Figure  6.2. 

In  all  three  starting  structures  it  was  observed  that 
most  of  the  spherodized  cementite  was  situated  at  the 
ferrite  grain  edges  (see  Figure  6.2).  The  volume  frac- 
tion of  the  cementite  particles  was  determined  for  the 
particles  situated  at  the  ferrite  grain  edges,  V^,  grain 
boundaries,  Vy,  and  inside  the  ferrite  matrix,  Vy,  in  the 
as-  coarsened  sample  of  batch  D^.  The  results  are  as 
follows : 

Vy  = 0.108 

Vy  = 0.03 

Vy  = 0.002 

Thus,  80%  of  the  spherodized  cementite  volume  was  situated 
at  the  ferrite  grain  edges.  It  may  be  concluded  that  the 
spherodized  cementite  is  situated  at  the  ferrite  grain 
edges  in  these  samples. 
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Figure  6. 


A photomicrograph  of  the  starting  structure 
D3  etched  with  the  double  etching  procedure. 
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A detailed  quantitative  microstructural  analysis 
was  carried  out  on  the  as-coarsened  structure  of  the 
batch  D^.  The  results  are  given  in  Table  6.17. 

Table  6.17 

Detailed  Microstructural  Analysis  of 
Starting  Structure 

1)  Volume  fraction  of  cementite  = 0.14 

2)  Surface  area  of  cementite  = 1384  cm  /cm'^ 

3)  Total  curvature  of  cementite  = 3.94  x 10°  cm/cm° 

4)  Surface  area  of  the  ferrite  grain  boundaries  = 

1744  cm2 / cm3 

5)  Length  of  the  ferrite  grain  edges  = 

2.024  X 106  cm2 / cm3 

6)  Length  of  the  ferrite  grain  edges  occupied  by 
cementite  = 2.26  x 10^  cm/cm3 

7)  Fraction  of  the  ferrite  grain  edge  length  occupied 

by  cementite:  0.53 

8)  Length  of  ferrite-ferrite-cementite  triple  lines  = 
4.171  X 10^  cm/ cm3 

It  is  observed  that  although  most  of  the  cementite  parti- 
cles are  on  the  ferrite  grain  edges,  only  537o  of  the 
ferrite  grain  edge  length  is  occupied  by  the  spherodized 
cementite.  Some  serial  sectioning  was  done  on  the  as- 
coarsened  structure  of  the  batch  . It  was  observed 
that  the  cementite  particles,  which  looked  discrete  and 
separate  on  the  section,  were  interconnected  in  the  third 
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dimension.  No  quantitative  measurements  were  made  to 
evaluate  the  connectivity  of  the  structure. 

A detailed  quantitative  microstructural  analysis  was 
performed  on  each  sample  of  one  series  (starting  struc- 
ture D^)  partially  austenitized  at  750°C.  There  are 
three  phases  present  in  the  partially  austenitized 
samples:  ferrite,  cementite  and  austenite  (which  is 
martensite  in  the  quenched  samples) . Let  superscript 
a denote  the  ferrite  phase , y denote  the  austenite  and 
C denote  the  cementite.  The  following  geometrical  prop- 
erties were  measured  in  each  sample  reacted  at  750  C. 

1)  Volume  fraction  of  austenite 

n 

2)  Volume  fraction  of  cementite  V^ 

3)  Total  surface  area  of  austenite 

4)  Total  surface  area  of  cementite 

n 

5)  Total  curvature  of  cementite 

6)  Surface  area  of  austenite-cementite  interface 

7)  Surface  area  of  ferrite-cementite  interface  Sy~°‘ 

8)  Surface  area  of  ferrite  grain  boundaries  Sy“ 

9)  Length  of  the  ferrite  grain  edges  Ly°“^ 

10)  Length  of  ferrite- ferrite-austenite  triple  lines  l““^ 

11)  Length  of  ferrite-ferrite-carbide  triple  lines  Ly“*^ 

12)  Length  of  ferrite-austenite-carbide  triple  lines  Ly^*^ 
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To  the  best  knowledge  of  the  author  such  an  exhaustive 
quantitative  microscopic  study  of  a phase  transformation 
has  never  been  done  before.  It  will  be  shown  in  the 
subsequent  sections  that  these  measurements  are  essential 
for  the  rigorous  understanding  of  this  apparently  simple 
process . 

In  order  to  etch  all  the  three  phases  and  the  ferrite 
grain  boundaries  present  in  the  partially  austenitized 
samples,  the  samples  were  first  etched  by  using  the 
etching  procedure  described  in  Table  6.2  and  then  re- 
etched by  using  the  procedure  described  in  Table  6.16. 

The  typical  microstructures  of  the  partially  austenitized 
samples  are  shown  in  Figures  6.3  to  6.5.  The  detailed 
quantitative  microstructural  analysis  of  the  samples 
reacted  at  750°C  was  performed  manually.  The  measure- 
ments were  carried  out  by  projecting  the  image  on  the 
T.V.  screen  of  the  Quantimet.  A grid  was  superimposed  on 
the  image  to  perform  the  measurements.  The  results  are 
given  in  Tables  6.18  to  6.20.  The  following  qualitative  observa- 
tions were  made  on  the  partially  austenitized  samples. 

(1)  The  spherodized  cementite  was  not  completely 
covered  with  austenite  til  the  austenite  volume  fraction 
was  approximately  0.5. 

(2)  The  growth  of  austenite  was  inhibited  across  the 


ferrite  grain  boundaries. 
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Figure  6.3.  Microstructure  of  a sample  partially  austeni- 
tized at  750°C  for  35  sec  (starting  structure 
D3). 
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Figure  6.4.  Microstructure  of  a sample  partially  austeni- 
tized at  750°C  for  50  sec  (starting  structure 
D3)  . 
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Figure  6. 


5.  Microstructure  of  a sample  partially  austeni- 
tized at  750°C  for  90  sec  (starting  structure 
D3). 
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Table  6.18 

A Detailed  Stereological  Analysis  of  the 
Interfacial  Areas 
Austenitizing  temperature:  750°C 

Starting  structure: 

2 

Units  for  surface  areas:  cm  /cm 


Austenitizing 
t ime  (sec) 

Sy  total 

r 

Sy  total 

cY-C 

cY-ct 

cOt-C 

cOiCt 

0 

0 

1384 

0 

0 

1384 

1744 

30 

968 

1162 

344 

624 

816 

1105 

35 

1042 

1016 

300 

742 

716 

911 

40 

1317 

995 

400 

917 

596 

744 

45 

920 

980 

■k 

JU 

k 

779 

50 

1789 

850 

491 

1298 

359 

668 

65 

1950 

681 

550 

1400 

131 

320 

90 

1763 

630 

606 

1157 

0 

258 

150 

1051 

400 

400 

651 

0 

108 

* 


No  measurements  performed 
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Table  6.19 

A Detailed  Stereological  Analysis  of  the 
Lineal  Features 

Austenitizing  temperature:  750°C 

Starting  Structure; 

3 

Units : cm/ cm 


Austenitizing 
time  (sec) 

j aaa 

f aay 

■r  aaC 

T ctyC 

0 

2.024x10^ 

0 

4.171x10^ 

0 

30 

1.05x10^ 

9.58x10^ 

1.851x10^ 

7.42x10^ 

35 

0.91x10^ 

8.7x10^ 

1.085x10^ 

4.81x10^ 

40 

0.694x10^ 

9.0x10^ 

0.85x10^ 

6.05x10^ 

45 

0.72x10^ 

■k 

* 

4.67x10^ 

50 

0.688x10^ 

1.14x10^ 

0.8x10^ 

9.05x10^ 

65 

0.237x10^ 

1.0x10^ 

*>v 

Vc 

90 

0.188x10^ 

9.32x10^ 

0 

0 

150 

0.086x10^ 

4.15x10^ 

0 

0 

No  measurements  made 
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Table  6.20 

Measured  Volume  Fractions  of  Different  Phases 
Austenitizing  temperature:  750°C 

Starting  Structure: 


Austenitizing 
time  (sec) 

V 

a > 
> 

0 

0.14 

0 

0.86 

30 

0.11 

0.116 

0.77 

35 

0.12 

0.130 

0.75 

40 

0.113 

0.184 

0.70 

45 

0.10 

0.160 

0.74 

50 

0.0847 

0.285 

0.63 

65 

0.0612 

0.45 

0.49 

90 

0.0548 

0.58 

0.37 

150 

0.037 

0.83 

0.13 
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(3)  Austenite  nucleated  perf erentially  along  the 
ferrite  grain  edges. 

(4)  The  austenite  growth  along  the  ferrite  grain 
boundaries  was  much  more  than  that  inside  the  ferrite 
grains . 

(5)  There  is  some  incubation  time  involved  in  the 
formation  of  austenite;  the  magnitude  of  this  incubation 
time  depends  on  the  austenitizing  temperature.  Further- 
more, the  magnitude  oJ^the  incubation  time  is  higher  for 
X^e_coarje^r„starting  structures.  Some  ^2  samples  were 

reacted  at  730°C.  It  was  observed  that  the  austenite 
was  not  formed  in  these  samples  for  reaction  times  less 
than  75  seconds.  Thus,  it  can  be  concluded  that  the 
observation  of  incubation  time  is  not  an  artifact 
arising  from  the  experimental  procedures. 

6 . 4 Data  Processing: 

(1)  It  was  observed  that  the  voixjme  fraction,  V^,  of 

n 

cementite  varied  linearly  with  the  surface  area,  Sy,  of 

cementite.  This  behavior  is  shown  in  Figures  6.6  to  6.10 

for  different  austenitizing  temperatures  and  starting 

structures.  Table  6.21  gives  the  values  of  slope  and 

C r 

intercept  of  the  plot  of  Sy  vs . Vy  for  different 
austenitizing  temperatures  and  starting  structures.  It  is 
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Figure  6.6.  A plot  of  the  surface  area,  Sy,  of 
cementite  vs.  its  volume  fraction, 
during  the  austenitization  of 
the  starting  structure,  D]^ , at  750°C. 
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Figure  6.7.  A plot  of  the  surface  area,  Sy,  of 

cementite  vs.  its  volume  fraction,  Vy , 
during  the  austenitization  of  the 
starting  structure,  D2 , at  740°C. 
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Figure  6.8.  A plot  of  the  surface  area,  Sy,  of 

cementite  vs.  its  volume  fraction,  Vy , 
during  the  austenitization  of  the 
starting  structure,  D2 , at  750°C. 
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Figure  6.9.  A plot  of  the  surface  area,  Sy,  of 

cementite  vs.  its  volume  fraction,  Vy , 
during  the  austenitization  of  the 
starting  structure,  D2 , at  760°C. 
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Figure  6.10.  A plot  of  the  surface  area,  Sy , 

of  cementite  vs.  its  volume  frac- 
tion, , during  the  austenitiza- 
tion of  the  starting  structure, 

D3 , at  750°C. 
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Aus  tenitizing 
temperature 

Table 
= A 

by 

Starting 

structure 

6.21 
v5  + A2 

Slope  A^ 

2/  3 

cm  /cm 

Intercept 

2 

cm  / cm 

750°C 

‘'i 

1.7x10^ 

70 

750°C 

D2 

1.12x10^ 

40 

750°C 

0.95x10^ 

60 

740°C 

D2 

1.21x10^ 

40 

760°C 

D2 

1.12x10^ 

40 

interesting  to  note  that  the  slope  and  the  intercept  do 

not  vary  significantly  with  the  austenitizing  temperature 

but  depend  on  the  starting  structure.  This  iriay_  indicate 

_that  __the  linear  relationship  between  Sy  and  Vy  results 

from  the  geometry  of  the  starting  microstructure  rather 

than  the- kinetics _of  the  process.  Figure  6.11  shows  a 

plot  of  (Sy/Sy^)  vs.  (Vy/Vy  ) for  all  the  austenitizing 

times,  temperatures  and  the  starting  structures  of 

V r 

spherodized  cementite  (V..  and  St,  ) are  the  volume  fraction 
and  surface  area  of  spherodized  cementite  in  the  unreacted 
samples.  All  the  data  points  fall  on  one  straight  line. 
Thus  we  have  the  following  empirical  relationship  between 
the  volume  fraction  and  surface  area  of  cementite  during 
the  austenitization  of  these  spherodized  structures : 


o < 
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Figure  6.11.  A plot  of  (Sy/SvQ)  vs. 

(Vy/VyQ)  during  the  austeni- 
tization of  the  spherodized 
cementite  structures . 


(6.6) 


+ 0.03 


V 


V, 


0 


Equation  (6.6)  is  valid  for  all  the  austenitizing  tem- 
peratures and  starting  structures  of  spherodized 
cementite.  The  mean  linear  intercept,  , of  the 
cementite  phase  is  given  by  the  following  equation: 


= 


4V 


V 


’V 


Combining  equations  (6.6)  and  (6.7)  give: 


A^  = 


A ^ 
^0 


V 


^0 

[1.01+.03— 


(6.7) 


(6.8) 


where  Aq  is  the  mean  intercept  of  spherodized  cementite  in 
the  starting  structure.  It  is  interesting  to  note  that, 
for  > 0.03, 


^ (6.9) 

Thus  mean  intercept  of  the  spherodized  cementite  does  not 
vary  significantly  with  the  austenitization  time  til  the 
process  is  almost  complete.  This  is  a rather  peculiar 
result,  because  one  would  expect  the  mean  intercept,  A^ , 
to  decrease  systematically  with  time  as  the  spherodized 
cementite  dissolves. 
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n 

It  was  observed  that  the  total  curvature,  M^,  of  the 
spherodized  cementite  varied  linearly  with  its  surface 

Q 

area,  S^,  during  the  austenitization  process.  This  be- 
havior is  shown  in  Figures  6.12  to  6.16,  for  different 
austenitizing  temperatures  and  starting  structures  of 
spherodized  cementite.  Table  6.22  gives  the  values  of 
the  slope  and  intercept  of  the  plot  of  vs.  for 
different  austenitizing  temperatures  and  starting  struc- 
tures . 


Table  6.22 
= A 4-  A 


Austenitizing 

temperature 

Starting 

structure 

Slope 

1/cm 

Intercept  A^ 
cm/cm^ 

750 

C 

^1 

4900 

4.0x10^ 

750 

C 

D2 

3000 

0.50x10^ 

750 

C 

“3 

2790 

0.50x10^ 

740 

C 

D2 

3180 

-2.0x10^ 

760 

C 

°2 

3000 

0.5x10^ 

Figure  6 

.17  shows  a plot  of  ) 

V Vq 

vs . ( 

(sj/sf;  ) for  all 

V Vq 

the  austenitizing  times,  temperatures 

and 

the  starting 

r c 

structures  of  the  spherodized  cementite  (M„  and  S 

''o  Vo 

the  total  curvature  and  surface  area  of  the  cementite 


are 
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tization  of  the  starting  structure,  Di , 
at  750°C. 


Figure  6.12. 
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Figure  6.13.  A graph  of  My  vs.  Sy  during  the 
austenitization  of  the  starting 
structure  D2  at  740°C. 
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Figure  6.14.  A graph  of  My  vs.  Sy  during  the 
austenitization  of  the  starting 
structure  D2  at  750°C. 
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c c 

Figure  6.15.  A graph  of  My  vs.  Sy  during  the 
austenitization  of  the  starting 
structure  D2  at  760°C. 
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Figure  6.16.  A graph  of  My  vs.  Sy  during  the 
austenitization  of  the  starting 
structure  D3  at  750°C. 
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o 


A graph  of  (My/Mvp)  vs. 
(SS/S^q)  during  tne  austeniti- 
zation of  the  spherodized 
cementite  structures . 


Figure  6.17. 
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phase  in  the  starting  structures) . The  plot  is  linear 
and  all  the  data  points  fall  on  one  straight  line.  Thus 
we  have  the  following  empirical  relationship  between  the 
total  curvature  and  surface  area  of  the  spherodized 
cementite  during  the  process: 


M 


V 


M. 


^0 

The  mean  curvature , 
given  by  the  following 


= 0.98[ 


] + 0.02 


(6.10) 


, of  the  spherodized  cementite  is 
equation  [34,35] : 


6.11) 


Combining  equations  (6.10)  and  (6.11)  yields 

= H^[0.98  + 0.02-^]  (6.12) 

Sy 

where  Hq  is  the  mean  curvature  of  the  spherodized 
cementite  in  the  starting  structure.  Thus  the  mean  curva- 
ture of  the  spherodized  cementite  does  not  vary  signifi- 
cantly with  the  austenitizing  time  til  the  process  is 
almost  complete.  Again,  this  is  a very  peculiar  result 
because  one  would  expect  to  increase  systematically 
with  the  austenitizing  time. 

Define  an  inverse  cumulative  frequency  function. 
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F^(£,t)  for  the  linear  intercepts  of  spherodized  cementite 
at  time  t as 

N (£,t) 

L 

where  Nj^(t)  is  the  total  number  of  intercepts.  Note 
that 

sH(t) 

NL(t)  = — (6.14) 

Figures  6.18  to  6.22  show  the  plots  of  F^(£,t)  vs. 

different  austenitizing  times  for  given  austenitizing 
temperatures  and  starting  structures.  It  is  observed 
that  for  a given  austenitizing  temperature  and  starting 
structure,  the  plot  of  F^(£,t)  vs.  (£/X^)  is  independent 
of  the  austenitizing  time.  In  each  of  the  figures  (6.18 
to  6.22),  it  is  observed  that  the  data  points,  at  dif- 
ferent austenitizing  times,  fall  on  the  same  curve. 

Thus  , 


F>(£,t)  = fu^txon  ofj  ^ <i)(£/A^)  (6.15 

(I/a)  only 

This  is  a very  peculiar  observation,  because  if  it  is 
due  to  the  kinetics  of  the  process , then  it  suggests  an 
asymptotic  distribution  of  the  cementite  particles  during 
the  austenitizing  process,  which  is  not  at  all  expected. 
On  the  other  hand  this  observation  can  be  explained  in  a 
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Figure  6.18.  A graph  of  inverse  cumulative 
frequency  function,  F>(£),  vs. 
l/X,  at  various  times  during  the 
austenitization  of  the  starting 
structure  at  750°C. 
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Figure  6.19.  A graph  of  inverse  cumulative 
frequency  function,  F>(£),  vs. 
£/A,  at  various  times  during  the 
austenitization  of  the  starting 
structure  D2  at  740°C. 
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Figure  6.20.  A graph  of  inverse  cumulative  frequency- 

function,  F>(£),  vs.  £/A,  at  various  times 
during  the  austenitization  of  the  starting 
structure  D2  at  750°C. 
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Figure  6.21.  A graph  of  inverse  cumulative 
frequency  function,  F>(£),  vs. 
£/A.,  at  various  times  during  the 
austenitization  of  the  starting 
structure  D£  at  760°C. 
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Figure  6.22.  A graph  of  inverse  cumulative 
frequency  function,  F>(£),  vs. 
l/\,  at  various  times  during  the 
austenitization  of  the  starting 
structure  D3  at  750°C. 
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very  simple  manner  in  terms  of  the  geometry  of  the 
starting  structure. 


6.4.1  Data  Processing  of  the  Detailed  Quantitative 

Microscopic  Measurements  Ilade  on  the  D-:^  Samples 
Austenitized  at  750”U 

(1)  Time  variation  of  the  volume  fraction  of 
austenite : 

A plot  of  log{VJ/l-V^}  vs . log  t is  shown  in  Fig- 
ure 6.23.  The  plot  linear  and  the  straight  line  has  a 
slope  of  2.25.  Thus  the  following  equation  describes  the 
time  variation  of  the  volume  fraction  of  austenite: 


V. 


1-V 


V 


(6.16) 


where  b]^  = 5.75  x 10  ^ per  (sec)^'^^  at  750°C.  Note  that 
B^  is  expected  to  change  with  the  austenitizing  tempera- 
ture. The  conventional  plot  of  log{ln(l/ (1-V^) } vs.  log  t 
was  found  to  be  non-linear  which  may  be  due  to  the  non- 
random  impingement.  This  is  not  surprising,  because 
austenite  is  not  randomly  distributed  in  the  matrix  and 
it  grows  preferentially  along  the  ferrite  grain  bounda- 
ries . 


(2)  Relationship  between  the  volume  fraction  and 
surface  area  of  austenite: 

A plot  of  {Sy/(1-Vy)}  vs.  is  shown  in  Figure  6.24. 
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Logt 

Figure  6.23.  A graph  of  Log(Vv/l-Vv) 
vs.  Log  t during  the 
austenitization  of  the 
starting  structure  D-?  at 
750°C. 
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zation  of  the  starting 
structure  D3  at  750°C. 
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The  behavior  is  linear.  The  following  equation  describes 
the  relationship  between  the  volume  fraction  of  austenite 
and  its  surface  area  during  the  process. 


V 

(1-V^) 


= B2  Vj  + B 


(6.17) 


where  B2  = 7100  per-cm  and  B^  = 250  per-cm.  Combining 
equations  (6.16)  and  (6.17),  together  with  the  definition 
of  the  surface  area  everaged  interface  velocity  of 
austenite,  Vg,  gives  the  following  result: 

.-4 


“I 


3.16x10 


cm/ sec 


(6.18) 


6 . 5 Data  Analysis: 

It  was  observed  that  80%  of  the  spherodized  cemen- 
tite  volume  was  situated  at  the  ferrite  grain  edges . 

This  spherodized  cementite  can  be  in  the  form  of 

(1)  complete  grain  edge  network, 

(2)  discrete  particles  situated  at  the  ferrite 
grain  edges 

or 

(3)  incomplete  grain  edge  network. 

Only  534  of  the  ferrite  grain  edge  length  was  occupied  by 
the  spherodized  cementite,  so  the  spherodized  cementite 
can  not  be  in  the  form  of  a complete  grain  edge  network. 
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Thus,  the  first  case  (i.e.  the  possibility  of  complete 
grain  edge  network)  can  be  rejected.  Let  us  consider  the 
second  possibility.  If  the  spherodized  cementite  is  in 
the  form  of  discrete  particles  of  the  same  shape,  then 
the  three  dimensional  particle  size  distribution  can  be 
evaluated  as  follows. 

Combining  equations  (6.13),  (6.14)  and  (6.15)  gives 


s^(t) 


(})  [5,/x^; 


By  definition 


nj^(5,,t)  = - 


dN^^(£,t) 


Thus 


(6.19) 


n.  (£ , t)  = 


s$(t) 


1 d{4)(&/A^) 
d{£/I^} 


Combining  equations  (6.20)  and  (2.10)  yields 


riy  (^  . t) 


T * 2k 


(6.20) 


{A^(t)} 


d(£/A^)  ' ^ 


} (6.21) 


d(£/A^) 


where  n^(£,t)  is  the  P.S.D.  function  and  K is  a shape 
factor.  The  second  factor  on  the  right-hand  side  of  this 
equation  is  a function  of  (£/A^)  only.  Let 


0 


[£/A^]  = ^ 

d(£/A^)  ^ . r. 


d(£/A^) 


(6.22) 


Then 
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n^(£,t)  = 


^ • 0(£/X^) 


(6.23) 


Let  us  define 


0>(«'/A^)  = / Q(l/X^)  • d(£/A^)  (6.24) 

l/X^ 

Integration  of  equation  (6.23)  gives  the  following 

expression  for  the  total  number  of  cementite  particles 

r 

per  unit  volume,  N^(t) : 


C 0^(0)  • S^(t) 

N^(t)  = ^ 


00 

where  0 (0)  = / Q(l/X^)  • d(£/X*^)  . Combining  equa- 

0 

tions  (6.25),  (6.8)  and  (6.6)  yields 

X 


(6.25) 


r 0>(O)Sv„  \ 

Ny(t)  - Zc  2 


''0  2 ''v  3 

C [-^]^[1.01(-^)+0.03]-^ 

XT'-" 

\ 


(6.26) 


V c 

During  most  of  the  transformation  V^/V^  >>  0.03.  Thus 


one  can  approximate 


. 40  (O)V^(t) 

N^(t)  - . ^ ^ 


71^ 


(6.27) 


For  a given  starting  structure,  Xg  is  constant  because 
it  is  the  initial  value  of  the  mean  intercept . Thus  equa- 
tion (6.27)  predicts  that  the  total  number  of  cementite 
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particles  per  unit  volume  is  directly  proportional  to 
its  volume  fraction  during  the  austenitization  process. 
Equations  (6.6),  (6.10)  and  (6.27)  give  the  following 
result . 


C 40  ro) 

N^(t)  = 


\(t)  = { 


40>(O)  \ 

2-H(^)s5(t)  - 0.355$  } 
^0 


yC  gC 

40^(0)  \ S 

= c 3 >H-^)(-7^)M$(t)  - .055$  } 

f^0>  ''  V, 


(6.28) 


5ince  A«,  0 (0) , 5$  , V$  and  M$  are  constants  for  a 
U > Vq  Vq  Vq 

given  starting  structure,  equation  (6.28)  predicts  linear 
relationships  between  N$(t) , V$(t) , 5$(t)  and  M$(t) . For 
a system  of  discrete  particles,  this  would  imply  the 
following . 

The  growth  rate  of  austenite  is  extremely  large,  so 
that  as  soon  as  the  austenite  grain  nucleates  on  a given 
cementite  particle,  the  cementite  particle  dissolves  com- 
pletely in  a very  short  time.  Thus  the  process  of 
cementite  dissolution  simply  consists  of  removing  the 
cementite  particles  from  the  system  at  a finite  rate.  If 
this  is  so,  then  5$  ^ should  be  a very  small  fraction  of 
the  total  cementite  surface  area  at  any  given  time.  This 
contradicts  with  the  experimental  observation  (see 
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Table  6.18)  that  for  > 0.28,  S^~*"/Sy  > 0.5.  Thus  the 
assumption  that  the  spherodized  cementite  is  in  the  form 
of  discrete  particles  of  the  same  shape  (at  any  given 
time)  leads  to  a contradiction. 

If  it  is  assumed  that  the  spherodized  cementite  is  in 

a form  of  discrete  particles  but  all  the  particles  do  not 

have  the  same  shape  at  any  given  time,  then  it  is  not 

possible  to  evaluate  the  three  dimensional  particle  size 

distribution  from  the  distribution  of  intercepts  and 

equation  (6.27)  is  not  valid.  However,  linear  relation- 
C C C 

ships  between  V^,  Sy  and  are  an  experimental  observa- 
tion and  this  behavior  cannot  be  explained  if  the  parti- 
cles dissolve  by  a three  dimensional  dissolution  geo- 
metry. It  may  be  possible  to  explain  a part  of  the 
experimental  behavior,  it  if  is  further  assumed  that  the 
cementite  particles  have  a long,  thin  shape  and  they  dis- 
solve in  a one  dimensional  fashion.  However,  this 
assumption  is  not  much  different  from  the  assumption  that 
cementite  is  in  the  form  of  an  incomplete  grain  edge 
network.  Furthermore,  if  it  is  assumed  that  "spherodized" 
cementite  is  in  the  form  of  an  incomplete  grain  edge  net- 
work, then  it  is  possible  to  explain  all  the  experimental 
observations  in  a simple  manner.  Thus  it  is  reasonable 
to  postulate  that  the  "spherodized"  cementite  is  present 
in  the  form  of  an  incomplete  grain  edge  network. 
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Since  the  network  is  along  the  ferrite  grain  edges , 

the  initial  length  of  ferrite-ferrite-cementite  triple 
• cx  cx  C 

lines,  L„  , should  be  equal  to  three  times  the  length  of 
° C 

the  network,  . Thus, 

''o 


, aaC  _ ^-r  C 


(6.29) 


The  initial  total  curvature  of  such  a network  is  given 
by  [122] 


M = ttL^ 


(6.30) 


Combining  equations  (6.29)  and  (6.30)  gives 


M 


Vr 


■ aaC 

"Vo 


(6.31) 


For  the  initial  structure  of  batch  (see  Table  6.16) 

mH  = 3.94  X 10^  cm/ cm^ 

and 


jCtotC  1 in6  / 3 

Ltt  - 4.17  X 10  cm/ cm 

'^O 

Thus  equation  (6.31)  is  satisfied. 

The  geometrical  description  of  cementite  as  an  in- 
complete grain  edge  network  explains  observations 
made  on  the  starting  structures.  These  observations  are 
as  follows . 

(1)  Almost  all  the  cementite  particles  are  observed 
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on  the  ferrite  grain  edges,  and  quantitatively  30%  of  the 
cementite  volume  is  situated  along  the  ferrite  grain 
edges . 

(2)  Although  the  cementite  is  situated  on  the 
ferrite  grain  edges,  only  53%  of  the  ferrite  grain  edge 
length  is  occupied  by  cementite. 

(3)  Cementite  particles  which  appear  discrete  on  a 
section  are  interconnected  in  the  third  dimension  (quali- 
tative serial  sectioning  evidence) . 

“Vg  “ 1-5^ 

This  geometrical  description  will  now  be  used  to 
explain  the  observations  made  on  the  cementite  phase 
during  the  austenitization  process. 

Since  the  cementite  is  assumed  to  be  in  the  form  of 
an  incomplete  network,  its  maximum  dimension  along  a 
cross-section  may  be  assumed  to  be  much  less  than  the 
length  of  each  separate  part.  (Note  that  this  is  not  an 
independent  assumption.)  It  will  be  further  assumed  that 
the  cementite  network  dissolves  in  a one  dimensional 
manner , such  that  its  length  decreases  with  time  while 
the  cross-sectional  area  remains  essentially  constant. 

Let  L^(t)  be  the  length  of  the  network  at  time  t.  The 
volume  fraction  of  cementite,  V^(t) , is  given  by 

C 

Vy(t)  = {Average  cross-sectional  area}  • L^(t)(6.32) 
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During  most  of  the  transformation,  the  surface  area, 

Q 

S^(t) , is  given  by 

Q 

Sy(t)  = {Average  cross-sectional  perimeter}L^(t)  (6. 33) 


During  a one  dimensional  dissolution,  the  average  cross- 
sectional  area  and  the  average  cross-sectional  perimeter 
does  not  change  with  time.  Thus, 


V. 


V. 


V 


0 


(Average  cross-sectional  area) 
(Average  cross-sectional  perimeter) 


(6.34) 


Combining  equations  (6.32)  to  (6.34)  yields 

S^(t)  _ V^(t) 

(6.35) 

Equation  (6.35)  predicted  by  the  model  is  very  similar  to 
the  empirical  equation  (6.6)  which  describes  the  experi- 
mental behavior.  The  only  difference  between  the  equa- 
tion (6.35)  and  (6.6)  is  the  non-zero  intercept  of  0.03. 
This  is  not  surprising  because  equation  (6.35)  is  not 
valid  near  the  end  of  the  transformation. 

DeHoff  [32]  has  shown  that  the  total  curvature  of  a 
tubule  network  is  directly  proportional  to  its  length. 

Thus  , 

M^(t)  - . Ly(t) 

where  is  a shape  factor. 


(6.36) 


182 


From  equations  (6.33)  and  (6.36) 


0 _ {Average  cross-sectional  oerimeter} 

V * (0.3/) 

0 

Combining  equations  (6.33),  (6.36)  and  (6.37)  yields 

M^(t)  s5(t) 

- = -V-  (6.38) 


M. 


V 


0 


0 


Equation  (6.38),  predicted  by  the  model,  is  very  similar 
to  the  empirical  equation  (6.10)  which  describes  the 
experimental  behavior.  The  only  difference  between  the 
two  equations  is  the  non-zero  intercept  of  0.02  in  the 
equation  (6.10).  Again,  this  is  not  surprising  because 
equation  (6.38)  is  not  valid  near  the  end  of  the  trans- 
formation . 

The  behavior  of  the  distribution  of  linear  inter- 
cepts of  the  cementite  phase  will  now  be  analyzed.  Since 
the  cementite  is  in  the  form  of  a tubule  network,  a large 
number  of  the  intersections  of  such  a network  with  the 
plane  of  polish  are  expected  to  be  along  its  cross- 
section  (i.e.  perpendicular  to  the  tubule  axis).  Thus, 
the  distribution  of  linear  intercepts  has  a major  contri- 
bution from  the  intercepts  from  the  cross-section  of  the 
network.  The  change  in  the  intercept  distribution  func- 
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P 


tion  with  the  austenitizing  time  is  due  to  the  change  in 
the  average  dimension  of  the  cross-section  and  the  change 
in  the  length  of  the  network  itself.  If  the  length  of 
the  network  is  decreased  simply  by  removing  some  parts  of 
the  network  (as  in  the  case  of  one  dimensional  dissolu- 
tion) , then  this  will  decrease  the  total  number  of  inter- 
cepts (the  decrease  will  be  proportional  to  decrease  in 
the  surface  area  of  the  network).  However,  as  the  length 
of  the  network  decreases , the  frequency  of  intercepts 
will  not  change  significantly.  Thus,  the  frequency 
function  is  expected  to  be  time  invariant  for  a given 
starting  structure.  The  observed  experimental  behavior 
is  thus  explained  by  the  model. 

6.5.1  Analysis  of  the  Global  Properties  of  Austenite 

It  is  observed  that  the  austenite  growth  is  inhibited 
across  ferrite  grain  boundaries.  This  may  be  due  to  an 
orientation  relationship  between  the  austenite  and  ferrite 
grains.  Judd  and  Paxton  [120]  found  that  the  austenite 
grains  have  {110}^  habit  plane.  The  same  orientation 
relationship  was  observed  by  Grozier,  Paxton  and 
Mullins  [123]  in  the  Fe-N  system.  The  presence  of  an 
orientation  relationship  between  the  austenite  and  ferrite 
grains  may  also  explain  the  observation  that  cementite  is 
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not  completely  covered  with  austenite  til  about  50% 
transformation.  It  was  also  observed  that  the  austenite 
grains  grow  preferentially  along  the  ferrite  grain 
boundaries.  These  observations  suggest  a plate-like 
growth  of  austenite  along  the  ferrite  grain  boundaries. 
Thus,  it  is  reasonable  to  assume  that  the  austenite  grows 
in  a two  dimensional  manner  along  the  ferrite  grain 
boundaries.  The  extended  volume  fraction  of  austenite, 
for  such  a two  dimensional  plate-like  growth,  is  given 
by 


where 


Vl  (t)  = A / A(t,T)N(T)dT 
ex  0 


(6.39) 


A = thickness  of  the  austenite  plates 
A(t,i)  = area  of  the  flat  surface  of  a plate 
nucleated  at  time  t 

t = total  time  for  which  the  transformation 
has  progressed 

N(t)  = extended  nucleation  rate  of  austenite 
Note  that  for  two  dimensional  growth  the  plate  thickness, 
A,  does  not  change  with  time.  Furthermore,  for  a given 
plate,  A is  expected  to  be  much  smaller  than  the  average 
dimension  along  the  flat  surfaces,  soon  after  nucleation. 
Thus , one  can  write  the  following  expression  for  the 
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extended  surface  area,  (t) , of  austenite: 

^ ex 

t 

si  (t)  = 2 / A(t,T)N(x)dT  (6.40) 

'^ex  0 

Combining  equations  (6.39)  and  (6.40)  gives 

(t)  = ^ S;^  (t)  (6.41) 

ex  ex 

The  relationship  between  the  real  volume  fraction,  V^(t) , 
and  surface  area,  Sy(t) , is  given  by  the  empirical  equa- 
tion (6.17)  and  it  is  as  follows: 


Sj(t)  = [l-Vj(t)] [B2V^(t)  + B3]  (6.15) 

where  B2  = 7100  cm^/cm^  and  B^  = 250  cm^/cm^.  Equa- 
tions (6.41)  and  (6.15)  can  now  be  used  to  evaluate  the 
impingement  function,  4^(Vy).  The  procedure  is  described 
in  section  3.5.  The  result  is  as  follows: 

lU(vl)  = (6.42) 

Also , 


Thus , 


dvj  - dV> 

ex 


l-vj(t) 


Vj  (t) 
ex 


(6.43) 


S^(t)  = (t) 

ex 


Furthermore , 
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Hence , 


(6.44) 


Using  equations  (6.43)  and  (6.44)  together  with  the 
definition  of  the  surface  area  averaged  interface  velocity 
gives 

In  order  to  make  quantitative  predictions  regarding  the 
time  variation  of  the  global  properties  of  austenite,  a 
kinetic  model  will  now  be  developed.  The  assumptions  are 
as  follows. 

(1)  The  nucleation  rate  of  austenite,  N(t) , can 
be  described  by  the  following  functional  form: 


where  b^^  is  the  kinetic  parameter  and  it  depends  on  tem- 
perature. The  exponent  m is  a constant.  For  a constant 
nucleation  rate,  m = 0.  If  the  nucleation  rate  increases 
with  time,  then  m > 0.  If  the  nucleation  rate  decreases 
with  time,  then  m < 0.  A large  negative  value  of  m would 
indicate  site  saturation. 

(2)  The  two  dimensional  growth  of  austenite  plates 


N(t)  = bj^t"^ 


(6.46) 


occurs  (by  definition)  due  to  increase  in  the  area  of 
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their  flat  surfaces.  In  order  to  have  a simple  but 
flexible  model,  the  following  functional  form  is  assumed 
to  describe  the  rate  of  change  of  the  area  of  flat  sur- 
faces, A(t,x),  of  a plate  nucleated  at  time  t. 

nb. 


‘•dt  ^ T 


G 

7^ 


(6.47) 


or 


A(t,x)  = A(t-x)  = bg/^(t-x)^/’^  (6.48) 

For  a plate  of  circular  shape  and  n = 1,  equation  (6.47) 
describes  semi-infinite  diffusion  controlled  growth.  Note 
that  bg  is  a kinetic  parameter  that  depends  on  tempera- 
ture. The  constant  n is  determined  by  the  rate  con- 
trolling mechanism. 

Combining  equations  (6.48),  (6.39)  and  (6.43)  gives 


Vy(t) 

(i-v;j(t)) 


{Ab^bp^^’^a 
N G mn 


}t 


(m+^l) 


(6.49) 


where 


= /u'"[l-u]^/"'du  (6.50) 

mn  Q 

Note  that  is  just  a number  and  its  value  depends  on  m 

and  n.  Let 


B 


Y 

1 


1 1/n 


mn 


(6.51) 
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Then, 


Vj(t)  _ V 


(6.52) 


Note  that  b|  is  independent  of  time  but  it  depends  on 
the  temperature. 

Equation  (6.52),  derived  from  the  model,  is  very 
similar  to  the  empirical  equation  (6.16)  which  describes 
the  experimental  behavior  (see  Figure  6.23).  For  compari- 
son, equation  (6.16)  is  rewritten  as  follows: 


where  b|  = 
equations 


i-vj(t) 


5.75  X 10'^  1/sec)^-^^  at  750°C. 
(6.52)  and  (6.16)  gives 


(6.16) 


Comparing 


m+i.  1.25 


(6.53) 


and,  at  750°C, 

“ 5.75  X 10'5  l/(sec)^'“  (6.54) 

Combining  equations  (6.48),  (6.40),  (6.44)  and  (6.53) 
gives 


s^rt) 

= qY  /.N  = /OK  ul/n  n^2.25 

-r  S„  (t)  (2bj^bg 


[l-Vl(t)]^  "^ex 


(6.55) 


Thus,  the  model  predicts  that  a plot  of  log  (t)  vs. 


ex 
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log  t should  be  linear  with  a slope  of  2.25.  Such  a plot 
for  the  experimental  data  is  shown  in  Figure  6.25.  The 
behavior  is  linear.  The  slope  of  the  straight  line  is 
2.16  which  is  very  near  to  the  value  of  2.25  predicted  by 
the  model.  The  intercept  of  this  plot  gives 

* 0-^2  (^.56) 

Substituting  equation  (6.56)  into  (6.54)  yields 

A=1.6xl0^cm  (6.57) 

It  is  now  possible  to  evaluate  the  surface  area  averaged 
interface  velocity  from  the  model.  Using  equa- 
tions (6.55),  (6.16)  and  (6.45)  and  the  definition  of  \lg 
gives 

77_1.8x10^  ! 

Vg  - ^ cm/sec  (6.58) 

Equation  (6.58)  is  very  similar  to  equation  (6.18)  ob- 
tained from  the  empirical  equations  (6.16)  and  (6.17). 

The  extended  length  of  the  aay  triple  lines,  (t) , 

''ex 

is  given  by 

L“''^(t)  = K / {A(t,T)}^/^N(x)dT  (6.59) 

ex  0 

where  is  a shape  factor,  determined  by  the  equation  of 
the  perimeter  of  the  flat  surfaces. 

If  the  perimeter  of  the  flat  surfaces  can  be  approxi- 
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Figure  6.25.  A graph  of  Log  vs. 

Log  t during  the  austeni- 
tization of  the  starting 
structure  D3  at  750°C. 
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mated  by  a circle,  then  is  equal  to  3.54;  for  a square 
perimeter  the  value  is  4.9.  Combining  equations  (6.59), 
(6.48)  and  (6.46)  gives 


aay 

"V 

ex 


(t)  = 


“m2n^^ 


where 


“m2n  ' /u'“[l-u)^''^"du 
Note  that  ajj^2n  J^st  a number.  Again, 

ex 

Thus  , 


(6.60) 


(6.61) 


(6.62) 


[l-Vj(t)]^ 


L7^(t) 

ex 


, , l/2n  , (“+2H+1) 

'Vn^G  m2n*'-  (6.63) 


The  model  predicts  that  a plot  of  Log{L““^(t)}  vs.  log  t 

ex 

should  be  linear.  Such  a plot  for  the  experimental  data 
is  shown  in  Figure  6.26.  The  plot  is  linear,  as  predicted 
by  the  model.  The  slope  of  the  plot  is  1.75.  Thus, 

m + ^ + 1 = 1.75  (6.64) 

The  intercept  of  the  plot  gives 

Equations  (6.64)  and  (6.53)  give  the  following  results. 
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Figure  6.26.  A graph  of  Log  (LVe^)  • 
Log  t during  the  austeni- 
tization of  the  starting 
structure  D3  at  750°C. 
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n = 1.0 
m = 0.25 

Since  the  values  of  m and  n are  now  known, 

“mn  = “0.25,1  = = 0.36 

^ 1 > 
“m2n  " “0.25,2  " = 0.5 


(6.66) 


(6.67) 


It  is  now  possible  to  evaluate  and  b^  if  the  value  of 
the  shape  factor,  Kj^,  is  known.  As  indicated  earlier, 
is  equal  to  3.54  for  a circle  and  4.9  for  a square. 

For  an  ellipse,  the  value  of  varies  from  3.54  to  7.95 
as  the  axial  ratio  varies  from  one  to  ten.  Thus,  the 
choice  of  the  value  of  does  not  seem  to  have  a signi- 

ficant effect  on  the  calculated  values  of  b^  and  bg. 
Accordingly,  it  will  be  assumed  that  is  equal  to  3.54. 
Using  this  value  of  together  with  the  equations  (6.65), 
(6.54)  and  (6.67)  gives 

bj^  = 1.8  X 10^ 

{ } (6.68) 

bg  = 5.7  X 10'^ 

All  the  parameters  of  the  model  are  now  known  and  their 
values  at  750°C  are  as  follows: 
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m = 0.25 


n = 1.0 


A = 1.6  X 10  ^ cm 


{ 


} 


(6.69) 


K. 


L 


3.54 


= 1.8  X 10^  (cm)  ^ • (sec) 
bg  = 5.7  X 10  ^ cm^/sec 


It  is  now  possible  to  predict  the  coverage  of  the  ferrite 
grain  boundaries  by  the  austenite  grains.  Since  the 
austenite  growth  occurs  in  a two  dimensional  fashion  along 
the  ferrite  grain  boundaries , the  "extended"  coverage  of 
the  ferrite  grain  boundaries  should  be  approximately 
equal  to  half  of  the  extended  surface  area  of  the 
austenite  plates.  Thus, 


Substituting  the  values  of  the  parameters  b.-,  b„  and 


2.25 


(6.70) 


a 


0.25,1 


2.25 


(6.71) 


or 


gua  _ gOta 


log{ 9 } = 2.25  log  t - 0.44 

(i-v;j)^ 


(6.72) 
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Logt 

a a aa 

A graph  of  Log  [SVQ-SV/ 
(l-V^)-^jvs.  Log  t during 
the  austenitization  of 
the  starting  structure  D3 
at  750°C. 


Figure  6.27. 


196 


*■ 


V 

cm~1 


Figure  6.28.  A graph  of  Ly  vs.  Sv  during 
the  austenitization  of  the 
starting  structure  D3  at  750°C. 
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The  model  predicts  that  a plot  of  log{S?,°‘-S?,°‘/ (1-vX)^}  vs. 

Vq  V V 

log  t should  be  linear  with  a slope  of  2.25  and  the  inter- 
cept should  be  -0.44.  Such  a plot  for  the  experimental 
data  is  shown  in  Figure  6.27.  The  plot  is  linear  with  the 
slope  of  2.26  and  the  intercept  of  -0.48.  Thus  the 
agreement  between  the  prediction  of  the  model  and  the  ex- 
perimental behavior  is  very  good. 

The  present  model  assumes  that  the  austenite  growth 
occurs  predominantly  along  the  ferrite  grain  boundaries. 

In  such  a case,  the  length  of  the  ferrite  grain  edges  con- 
sumed by  the  austenite  grains  should  be  directly  pro- 
portional to  the  area  of  the  ferrite  grain  boundaries 
consumed  by  austenite.  Thus,  the  length  of  the  ferrite 
grain  edges,  Ly  (t) , should  vary  linearly  with  the  surface 
area  of  the  ferrite  grain  boundaries,  Sy°‘(t),  during  the 
process.  Figure  6.28  shows  a plot  of  Ly““(t)  vs.  Sy'^(t) 
for  the  experimental  data.  The  behavior  is  linear,  as 
predicted  by  the  model. 

The  model  developed  in  this  section  explains  the 
experimentally  observed  behavior  during  the  austenitiza- 
tion process. 


6 . 6 Conclusions : 

(1)  "Spherodized"  cementite  is  in  the  form  of  an 
incomplete  grain  edge  network. 
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(2)  During  the  austenitization  process  the  dissolu- 
tion of  the  cementite  network  takes  place  in  a one 
dimensional  manner. 

(3)  Austenite  growth  occurs  along  the  ferrite  grain 
boundaries . 

(4)  Austenite  grains  grow  in  a two  dimensional 
fashion.  The  rate  of  growth  of  an  austenite  plate  can  be 
described  by  the  rate  of  change  of  the  area  of  its  flat 
surfaces.  The  mathematical  form  is  as  follows: 


/dA\  _ 1 

where  A is  the  area  of  the  flat  surface  of  an  austenite 

plate.  The  parameter  b^  is  independent  of  time  but  it 

depends  on  the  austenitizing  temperature.  At  750°C,  b^ 

-7  2 

IS  equal  to  5 . 7 x 10  cm  /sec. 

(5)  The  nucleation  rate  of  austenite  increases  with 
time  and  it  can  be  described  by  the  following  functional 
form : 


N(t)  = b^t^/^ 

where  the  kinetic  parameter,  b^,  is  a function  of  tempera- 
ture but  independent  of  time.  At  750°C,  b^,  is  equal  to 
1.8  X 10^  per  cm^  (sec)^'^^. 

(6)  The  approach  developed  in  the  third  chapter  has 
been  successfully  used  to  interpret  the  experimental  data. 
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It  is  interesting  to  note  that  the  approach  successfully 
explains  the  experimentally  observed  behavior  of  nine 
independent  s tereological  properties  in  a quantitative 
manner.  To  the  best  knowledge  of  the  author,  this  has 
never  been  done  before  and  it  clearly  demonstrates  the 
potential  and  practical  feasibility  of  the  approach. 
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CHAPTER  VII 

CONCLUSIONS  AND  SUMMARY 

(1)  It  is  possible  to  evaluate  the  global  properties 
of  an  evolving  microstructure  from  a given  description  of 
the  local  variables.  The  approach  provides  a means  for 

a systematic  analysis  of  the  quantitative  microscopic 
data  on  the  kinetics  of  phase  transformations. 

(2)  A flexible  but  simple  model  has  been  developed 
to  evaluate  the  global  properties  of  a growing  phase 
during  a phase  transformation. 

(3)  A technique  has  been  developed  to  deduce  the 
nucleation  rate  and  growth  rate  from  the  time  variation 
of  the  global  properties. 

(4)  The  effects  of  various  growth  geometries  on 

the  kinematics  of  phase  transformations  have  been  explored 
in  detail.  It  is  shown  that  the  hyperbolic  time 
dependence  of  the  average  interface  velocities  (observed 
in  recrystallization  studies)  may  be  just  due  to  the  two 
or  one  dimensional  nature  of  the  growth  process. 

(5)  Impingement  equations  have  been  developed  for 
non-random  impingement.  A phenomenological  approach  has 
been  developed  to  evaluate  the  impingement  function  for 
two  dimensional  and  one  dimensional  growth  processes. 
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(6)  A generalized  continuity  equation  has  been  de- 
rived . 

(7)  An  algorithm  is  developed  for  the  prediction 
of  the  evolution  of  the  size  distribution  during  a dis- 
solution process. 

(8)  The  theoretical  analysis  is  applied  to  the  ex- 
perimental data  of  other  research  workers  on  the  kinetics 
of  the  recrystallization  process.  Every  qualitative  and 
quantitative  experimental  observation  is  explained  by 
using  the  present  approach.  It  is  concluded  that  the 
recrystallized  grain  size  (as  reported  by  the  line  inter- 
cept count)  will  be  independent  of  the  annealing  tempera- 
ture if  the  recrystallized  regions  grow  in  a two  dimen- 
sional or  one  dimensional  manner. 

(9)  The  kinetics  of  austenitization  of  a pearlitic 
structure  was  studied  at  three  temperatures.  The  volume 
fraction  and  surface  area  of  austenite  was  measured  at 
various  times  during  the  process.  The  nucleation  rate 
and  growth  rate  of  austenite  was  deduced  from  the  ob- 
served time  variation  of  the  volume  fraction  and  surface 
area  of  austenite.  It  is  concluded  that  the  process  takes 
place  by  a constant  nucleation  rate  and  constant  growth 
rate  of  austenite. 

(10)  A detailed  quantitative  microscopic  investi- 
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gation  was  carried  out  on  the  austenitization  of  three 
different  spherodized  cementite  structures.  The  experi- 
ments were  done  at  three  austenitizing  temperatures . 

(11)  The  "spherodized"  cementite  is  shown  to  be 
present  in  the  form  of  an  incomplete  grain  edge  network. 
During  the  austenitization  process  the  cementite  network 
dissolves  in  a one  dimensional  manner.  The  austenite 
grows  by  a two  dimensional  growth  along  the  ferrite 
grain  boundaries.  The  local  growth  rate  of  austenite  can 
be  described  by  a parabolic  time  law.  The  nucleation 
rate  of  austenite  increases  with  the  austenitizing  time 
and  it  is  directly  proportional  to  t'‘. 

(12)  The  present  approach  is  successfully  applied 
to  explain  the  experimentally  observed  behavior  of  nine 
independent  stereological  properties  during  the 
austenitization  of  the  "spherodized"  cementite  structures. 
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